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ABSTRACT 

We evaluate the effects of environment and stellar mass on galaxy properties at 0.85 < z < 1.20 using 
a 3.6/Ltm-selected spectroscopic sample of 797 cluster and field galaxies drawn from the GCLASS survey. 
We confirm that for galaxies with LogM,/MQ > 9.3 the well-known correlations between environment 
and properties such as star-forming fraction (tsp), SFR, SSFR, D„(4000), and color are already in 
place at 2: 1. We separate the effects of environment and stellar mass on galaxies by comparing the 
properties of star-forming and quiescent galaxies at fixed environment, and fixed stellar mass. The 
SSFR of star-forming galaxies at fixed environment is correlated with stellar mass; however, at fixed 
stellar mass it is independent of environment. The same trend exists for the D„(4000) measures of 
both the star-forming and quiescent galaxies and shows that their properties are determined primarily 
by their stellar mass, not by their environment. Instead, it appears that environment's primary role is 
to control the fraction of star forming galaxies. Using the spectra we identify candidate poststarburst 
galaxies and find that those with 9.3 < LogM^/M© < 10.7 are 3.1 ± 1.1 times more common in 
high-density regions compared to low-density regions. The clear association of poststarbursts with 
high-density regions as well as the lack of a correlation between the SSFRs and D„(4000)s of star- 
forming galaxies with their environment strongly suggests that at z ~ 1 the environmental-quenching 
timescale must be rapid. Lastly, we construct a simple quenching model which demonstrates that 
the lack of a correlation between the D„(4000) of quiescent galaxies and their environment results 
naturally if self quenching dominates over environmental quenching at z > 1, or if the evolution of 
the self-quenching rate mirrors the evolution of the environmental-quenching rate at 2; > 1, regardless 
of which dominates. 

Subject headings: galaxies: clusters, high-redshift, stellar content, star formation 
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1. INTRODUCTION 

It has been known for a long time that the properties 
of galaxies depend strongly on both their stellar mass 
and their larger-scale environment. In the last decade, 
the SDSS and other spectroscopic surveys have revo- 
lutionized this field, providing high-quality data that 
has quantified these correlatio ns with high preci sion 
(see e.g., the review bv .Blanton fc Moustaka"sll2009[ ). It 
is now clear that in the local universe practically any 
property that can be used to characterize a galaxy shows 
some correlation with both its environment and stellar 
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more frequently hav i ng elliptical/SO rriorpho loKies (e.g., 
iBamford etlD[200llNair fc Abrahamll2010D . 

There is also evidence for a covariance between 
steUar mass and environment, with more mas- 
sive galaxies residing in denser enviro nments (e.g., 
iKauffmann et al.l 12004 iBaldrv et all [2006) . although 
this i ssue has not yet reac h ed com plete consensus, see 
e.g., Ivon der Linden et al.l ()2010l ) for an alternative 
result. Either way, it is clear that any potential covari- 
ance between stellar mass and environment needs to be 
isolated in order to determine their relative influence on 
the galaxy population. 

Recently. iPeng et"al1 (|2010[ ) performed such an analy- 
sis and demonstrated that the star formation quenching 
effects of stellar mass and environment are independent 
of each other in the SDSS and zCOSMOS datasets up 
to z ^ 0.6. They used this separability to construct 
an empirically-based model which implies that "mass 
quenching" (which we will refer to in this paper as "self 
quenching") is the primary mechanism that ends star 
formation in massive galaxies (LogM^^/M© > 10.6) in all 
environments at all redshifts; whereas "environmental 
quenching" is increasingly important for lower-mass 
galaxies at decreasing redshift. 

Although we continue to acquire higher-quality mea- 
surements of the correlations between galaxy properties 
and their stellar mass and environment - and even 
have an empirical framework for inte rpreting how the y 
quench the overall galaxy population ()Peng et al.ll2010D . 
we still have only circumstantial evidence about the 
physical processes that are responsible for quenching. 
An obvious next step is to examine the evolution of 
galaxy properties as a function of stellar mass and 
environment over a wide range of epochs to determine 
the frequency of quenching and constrain the relevant 
timescale, an important clue for identifying the physical 
processes involved. 

In the last five years extensive spectrosc opic surveys 
of the 2 ~ 1 universe such as G OODS (lElbaz et al.l 
20071: iPopesso et al."2011'). DEEP2 dCooper et al."2006l. 



2007112 008. 2010b; Griitzb auch et al.l,2 011b). zCOSMOS 
(llovi^no et al.l I2010t IPeng et al.l 120101 : Cucciat i et all 



20101 ) . ROLES (iLi et all 1201 Ih . VVDS fiGucciati et al.l 
200a iScodeggio et al l I2009D. as well as the HiZELS 



narrow-band survey ( Sobral et al.ll201ll ) have begun to 
quantify the effects of stellar mass and environment on 
galaxy properties. Photometric surveys a re now doing 
the same at redshifts up to z ^- ^ 2-3 (e.g., lOuadri et al.l 
120071: iGniFzbauch et al.il2011al: lOuadri et al.M201lD . So 
far, this body of work has been somewhat controversial, 
with studies disagreeing on whether environment is an 
important factor in the quenching of star formation, 
and whether the star-formation-rate- (hereafter, SFR) 
density relation seen at z still exists at z ^ 1. 

Although the reason for the disagreements has not 
been dissected in detail, much of it may be semantic in 
nature, and not due to errors in data analysis or sample 
bias. This is because some authors have studied the 
color-density relation, others have looked at the specific- 



star-formation-rate- (hereafter, SSFR) density relation, 
while others yet have studied the actual SFR-density 
relation. Some studies correlate the mean color, mean 
SFRs, and mean SSFRs of particular subsamples of 
galaxies with environment. Others compare the change 
in the fraction of star-forming galaxies (hereafter isp) 
with environment, where star forming has been defined 
using any of those three criteria. At some level all of 
these correlations have been analyzed under the moniker 
of the "SFR-density relation", even though they are 
in fact quite different quantities. As we will show in 
this paper, the mean properties of galaxies and the 
behave quite differently as a function of environment 
and stellar mass, and these different behaviors, as well 
as some selection effects, may be at the root of many of 
the past discrepancies. 

Another major challenge for field galaxy surveys that 
has been largely ignored is that environment itself is 
difficult to define, and that the choice of parameteriza- 
tion almost certainly affects the interpretation of any 
analysis. Typically 3rd or 5th nearest neighbor distance 
is used as a proxy for galaxy environment; however, 
because it counts only a few galaxies this metric is 
noisy and only accounts for the density of galaxies on 
small scales. Using a larger number of nearest neighbor 
galaxies provides a better measurement of the larger 
scale environment; however, this comes at the cost of 
losing small scale information. Furthermore, increasing 
the number of nearest neighbors substantially leads to 
problems with the edges of survey fields and requires 
throwing out a large fraction of the data. A compre- 
hensive discussion of the issues ass ociated with definin g 
galaxy environment is presented in lCooper et al.l (|2005| l. 

Rich galaxy clusters offer an interesting alternative for 
measuring the effect of environment on galaxy evolution 
as they cover a wide range of environmental densities. In 
particular, their cores are the most extreme high-density 
environments at any epoch, and therefore leave little 
ambiguity about the measurement of local environment. 
Of course, the tradeoff with using clusters is that only a 
small percentage of galaxies in the local or high-redshift 
universe live in clusters, and therefore the environmental 
effects seen there may not be generalizable to the galaxy 
population at large. Still, they offer the possibility to 
define the upper limits of how strong environmental 
effects can be at any epoch. 

A handful of clusters at z ^ 1 have been studied 
with re gard to enyironmental and mass ev o lution 
effect s. JPatel et all (l2009l |2011D: iKovama et all (l2008l 
f20T0l) : iRosati et al.l (I2009D: iStrazzullo et al.l (12011 



Vulca ni et all (|2010D : iDemarco et al.l ()2010bf ) have all 
found that the correlation between SFR and SSFR with 
environment is still preserved in c luster s around z ^ 

1^ In additi on, iGobat et al.l ()2008D , and iRettura et al.l 

(|2010l 1201 ID have also found that the population of 
quiescent galaxies in clusters at z ~ 1 is more evolved 
than quiescent galaxies in the field. Interestingly, there 
is some evidence that at z > 1.4 even th e cores of clusters 
are becoming increa singly active (e . g.. [Tran et al] 120101: 
Hilton et all 120101: iHavashi et al.l 120101: iKuiper et al.l 
2010; Hat ch et al.l 201 ll). although this effe c t is not seen 
in al l clusters (e.g., iTanaka et al.l 120101 : I Bauer et al.l 
1201 If) . Whether this is truly the epoch where the 
expected change in the star formation properties of 
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cluster galaxies actually occurs, or is simply a selection 
effect is still difficult to tell as our current knowledge is 
based on just a handful of clusters. 

Here we present the most comprehensive study of the 
relative effects of stellar mass and environment on both 
cluster and field galaxies at z ~ 1 to date using data 
from the Gemini Cluster Astrophysics Spectroscopic 
Survey (GCLASS, Pis: Wilson + Yee). GCLASS is a 
spectroscopic survey of 10 rich clusters at 0.85 < z < 
1.34 using the CMOS instruments on Gemini-North and 
Gemini-South. The GCLASS cluster sample has been 
drawn from the Spitzer Adaptatio n of the Red-s e quence 
Cluster Survey (SdARCS. see IMuzzin et all l2009bl : 
IWilson et aLll2009l: iDemarco et al.ll2010aD . In this paper 
we use 9 of the clusters, and examine the properties 
of cluster galaxies as a function of stellar mass and 
environment for the first time using a large sample of 
clusters. 

One of the major improvements the GCLASS data 
offer over other spectroscopic studies is that the spectro- 
scopic targets have been selected based on their observed 
IRAC 3.6/im flux. This is different from other z ~ 1 
field galaxy surveys such as DEEP2, which is an R-band 
selected survey (R < 24.1), and zCOSMOS, which is an 
Lband selected survey (I < 22.5 for zCOSMOS-bright). 
At z ^ 1 these produce rest-frame U- and B-band 
selected samples, whereas the 3.6/im selection from 
GCLASS produces a rest-frame H-band selected sample 
and is therefore much closer to a stellar-mass-selected 
sample. 

This paper is laid out in the following manner. In § 2 
we define the cluster sample, § 3 discusses the reduction 
of the spectroscopic data, and in § 4 we explain how 
physical quantities such as stellar masses and SFRs are 
extracted from the data. Readers interested in the main 
results can skip those sections and refer to § 5 where 
we look at the correlation of galaxy properties for the 
sample of galaxies with stellar masses (LogM^/M©) > 
9.3 as a function of environment and stellar mass, and 
§ 6 where we look at the same properties but holding 
galaxy type, environment, and stellar mass fixed in 
order to untangle their interdependence. In § 7 we 
discuss the population of poststarburst galaxies as a 
function of environment and stellar mass, and in § 8 we 
present a simple model for the evolution of the quiescent 
galaxy population. We conclude with a summary in § 9. 
Throughout this paper we assume a ACDM cosmology 
with ilm = 0.3, f^A = 0.7, and Hq = 70 km s^^ Mpc"^ 

2. GCLASS CLUSTER SAMPLE 

The GCLASS cluster sample is selected from the 42 
deg^ SpARCS optical/IR cluster survey. SpARCS de- 
tects cluste rs using the clust e r red- sequence method de- 
veloped by iGladders fc Yei (l2000l . [2005) for the Red- 
seque nce Cluster Surveys (RCSl & RCS2, see lYee et al.) 
[2007t .Gilbank et al. 2011). Those surveys select clus- 
ters as overdensities of galaxies in an R - z' color space 
because the R - z' color spans the 4000A break fea- 
ture in early- type galaxies to z ~ 1. SpARCS selects 
clusters using a z' - 3.6/im color selection, which spans 
the break at z > 1. The z' data for the SpARCS sur- 
vey were obtained using the CFHT and CTIO 4m tele- 
scopes and the 3.6/im data used in SpARCS is taken 
from the 50 deg^ Spitzer Wide-area Extragalactic Sur- 



vey (SWIRE, iLonsdale et al] I2003D . An overview of 
the SpARCS observations, field layouts, an d photometry 
is prese nted in the compa nion papers by IMuzzin et al.l 
(l2009bD and IWilson eTall (|20O9). The SpARCS cluster 
detec tion algorithm is a slightly modified version of the 
Gla dders &: Yed ( 120001120051) algorithm and is discussed 
in IMuzzin et all (|2008[ ). The SpARCS cluster catalog 
will be presented in a future paper (A. Muzzin et al., in 
preparation) . 

The 10 clusters selected for foUowup as part of 
GCLASS were chosen in the following manner. Clusters 
with red-sequence photometric redshifts between 0.85 
< z < 1.25 were ranked in order of increasing detection 
significance. From this list we hand-selected a sample 
that attempted to include the richest clusters, but also 
to have an even distribution both in photometric red- 
shift and right ascension (R.A.). The resulting sample of 
nine clusters are evenly divided into three redshift bins, 
0.85 < z < 1.0, 1.0 < z < 1.1, 1.1 < z < 1.25, and are 
generally the three richest clusters in SpARCS in their 
respective redshift range. Lastly, we included SpARCS 
J003550-431224 at z = 1.34, a cluster for which we had 
already some observations, (|Wilson et al.l 120091 ) as the 
10th cluster to extend the redshift range of the sample. 
The high-redshift of SpARCS J003550-431224 means the 
spectra available for it are limited to Xrest < 4050A. This 
wavelength coverage is too short for the spectral prop- 
erties to be analyzed in the same manner as the other 
clusters, so in the current paper we exclude SpARCS 
J003550-431224 and focus on the remaining 9 clusters. 
All statistics of the spectroscopy quoted throughout the 
remainder of the paper exclude the spectra from SpARCS 
J003550-431224. 

We note that there are richer clusters in SpARCS that 
were not included as part of GCLASS because of the 
redshift and R.A. distribution requirements and there- 
fore, formally GCLASS is not a richness-limited sample 
of clusters at 0.85 < Zphot < 1-25. The GCLASS clusters 
comprise 10/24 of the richest clusters at Zphot > 0.85 in 
SpARCS and therefore in general, should be considered 
a fair sampling of IR-selected rich clusters within that 
redshift range. Color images of the GCLASS clusters in 
the gz'[3.6/xm] bands are shown in Figure 1. 

As Figure 1 shows, even rich clusters a.t z 1 have a 
diverse range of morphologies, from roughly spherically- 
symmetric and centrally-concentrated with a clear cen- 
tral galaxy, to asymmetric with filamentary-like struc- 
ture and no clear central galaxy. Figure 1 highlights the 
importance of using a sample of clusters when studying 
the effect of the cluster environment on galaxy evolution, 
so as to average over the peculiarities of the assembly his- 
tory of individual systems. 

Given the complexity of the galaxy distribution in 
the clusters, the centroid for 8 of the 9 clusters has 
been defined as the position of the brightest cluster 
galaxy (BCG), where the BCG is defined as the bright- 
est 3.6/im source that has a spectroscopic redshift con- 
sistent with being a cluster member. As Figure 1 il- 
lustrates, the location of the BCG for those 8 clusters 
is close to the center of the luminosity-weighted red- 
sequence galaxy distribution in the observed z' - 3.6/xm 
vs. 3.6/im color-magnitude space, and hence makes a log- 
ical choice of centroid. In the remaining cluster, SpARCS 
J105111+581803 at z = 1.035, the BCG is substantially 
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offset from the center of the luminosity-weighted red- 
sequence galaxy distribution. For this cluster, the cluster 
center has been defined as the centroid of the luminosity- 
weighted red-sequence galaxy distribution (a discussion 
of ho w this is computed is presented in iMuzzin et alJ 
MM) 

3. SPECTROSCOPIC DATA 

The goal of the GCLASS spectroscopic observations 
was to obtain spectra for ~ 50 members in each of the 
ten clusters. Based on the cluster color-magnitude dia- 
grams and richnesses, we estimated this would provide 
a spectroscopic completeness of > 50% for galaxies with 
Ms.e/xm < M*3.6^ni + 1-0 (see § 4.6 for a complete dis- 
cussion of the spectroscopic completeness). Geometric 
slit placement restrictions make obtaining 100% com- 
plete spectroscopic samples unrealistic, so the sampling 
in GCLASS was chosen as a trade off between complete- 
ness and efficiency. The completeness is sufficient to de- 
termine high-quality velocity dispersions for each clus- 
ter (see G. Wilson in preparation) and also allows us to 
compute reliable completeness corrections as a function 
of both stellar mass and clustercentric radius that can be 
used for defining complete samples of galaxies for galaxy 
evolution studies. Although the completeness level per 
cluster is not as high as some previou s surveys of indi- 
vidual clusters at sim ilar redshift fe.g.. lTran et al.l[2007l : 
iDemarco et all [2001 . GCLASS uses substantially fewer 
masks per cluster which in turn permits the study of a 
larger sample of clusters and cluster galaxies. 

Full details of the GCLASS mask design process, and 
data reduction will be presented in a future catalog re- 
lease paper, here we outline the most relevant steps of 
the observational setup and data reduction. 

3.1. Target Selection and Mask Design 

Spectroscopic targets were prioritized using three cri- 
teria, which in order of importance were: 1) cluster- 
centric radius 2) observed z' - 3.6^m color, 3) 3.6^m 
flux. Highest priority was given to bright 3.6/im sources 
in the cluster core with colors close to the observed red- 
sequence, lowest priority was given to faint 3.6/im sources 
on the cluster outskirts with z' - 3.6/im colors 1.0 mag- 
nitude redder or 1.5 magnitudes bluer than the cluster 
red-sequence^^. As an illustration of how the slit priori- 
ties were assigned, in Figure 2 we plot the observed z' - 
3.6/im vs. 3.6/im color-magnitude diagram for SpARCS 
J10511H-581803. 

In total there were 12 priority levels for targets (PI 
- P12). PI - PIG were the primary selection criteria 
optimized for obtaining the maximum number of clus- 
ter redshifts, and Pll and P12 were included as "filler" 
slits to maximize the number of slits per mask as well 
as provide a sample of galaxies without a specific color 
selection to test for possible selection biases. P1-P4 were 
the highest priority slits, targeting galaxies within 550 
kpc of the cluster center and z' - 3.6/im colors close to, 
or on the observed cluster red sequence (see Figure 2). 
P5 were sources within R < 550 kpc that met no spe- 
cific color criteria but with detections in the MIPS 24/im 

The precise color prioritization varies from cluster-to-cluster 
by a few tenths of a magnitude around these numbers. The color 
selection was broader for the higher-redshift clusters due to the 
bluer rest-frame probed by z' - 3.6/jm 



band. These sources were given higher priority than filler 
slits because of the possibility that they may be cluster 
galaxies with colors substantially different than the pri- 
mary cluster color priorities. P6 - PIG were identical to 
PI - P5, but for galaxies at R > 55G kpc. The filler slits, 
Pll and P12, were galaxies in any location with z' < 22 
mag, and z' > 22 mag, respectively. 

Slitmasks were designed using custom software created 
by us to produce output files which could also be read 
by the gmmps software provided by the Gemini observa- 
tory. Masks contained between 15 and 59 slits (excluding 
alignment stars), with a median of 42 slits per mask. Ta- 
ble 1 lists the clusters with the total number of masks and 
slits per cluster. 

The greatest concern with the selection of the spec- 
troscopic targets is that the use of a z' - 3.6/im color 
cut could introduce some bias into the galaxy sample. 
A pure 3.6/im fiux-limited sample would be preferable; 
however, due to the high-redshift of the clusters such a 
selection would be dominated by low-redshift foreground 
galaxies and is not an efficient use of telescope resources. 

As Figure 2 shows, the color priority selection was cho- 
sen to be quite broad so that it would span the properties 
of even the most extreme red and blue cluster members, 
and only exclude the most obvious very low/high redshift 
interlopers. It is also important to note that the color 
criteria only provided priority levels for the mask design 
algorithm, not hard color cuts. Many filler slits (Pll & 
PI 2) that did not meet any color criteria were placed 
and we used these to test the if the color prioritization 
provides an optimum sampling of cluster members. Filler 
slits resulted in a total of 274 redshifts across all clusters, 
15 of which were cluster members. This suggests that the 
color selection prioritized cluster members well (422/437 
of the cluster members met the color criteria), but also 
provided an efficient use of overall telescope time (only 
15/274 galaxies that did not meet the color criteria were 
members). 

3.2. Observations and Data Reduction 

The spectroscopic masks were observed with GMOS- 
N and GMOS-S using the R15G grating with a central 
wavelength of 75GGA with 1" wide slits which gives a 
resolution of 17A, equivalent to an R = 44G. For all 
masks we used 3" long microslits with the nod-and- 
shuffie option (NfcS, see l azebrook & Bland-Hawthor^ 
l200lt [Abraham et"all [2001 available with CMOS. For 
each cluster, approximately half the masks were observed 
in the " band shuffle" mode, and the other half were ob- 
served with the "micro shuffle" mode. The main advan- 
tage of the band shuffle mode is that the shuffled charge 
is stored on the top and bottom third of the chip while 
observing, not next to the slit and therefore the 3" slits 
can be packed directly beside each other in the cluster 
core where the density of galaxies is highest. The disad- 
vantage is that because of the need to store the shuffled 
charge, only a strip covering the central 1.7' of y-axis the 
chip can be used. In the micro shuffle mode, charge is 
stored beside the slit, which means slits must be at least 
3" away from each other; however, the full 5.5' FOV of 
CMOS can be used. We found that the combination of 
half band shuffle and half micro shuffle masks produced 
a distribution of slits that provided both a good sam- 
pling of the cluster core, but also fair coverage of cluster 
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TABLE 1 

GCLASS Cluster Sample 



Name 




R.A. 


Decl. 


Masks 


Total 


Total 


Cluster 


Reference 






J2000 


J2000 




Slits 


Redshifts 


Membersll] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SpARCS J003442-430752 


0.867 


00:34:42.086 


-43:07:53.360 


4 


159 


137 


45 (40) 




SpARCS J003645-441050 


0.869 


00:36:45.039 


-44:10:49.911 


4 


170 


119 


47 (47) 




SpARCS J161314+564930 


0.871 


16:13:14.641 


56:49:29.504 


5 


211 


161 


92 (85) 


Demarco et al. (2010a) 


SpARCS J104737+574137 


0.956 


10:47:37.463 


57:41:37.960 


4 


200 


147 


31 (31) 




SpARCS J021524-034331 


1.004 


02:15:23.200 


-03:43:34.482 


4 


167 


125 


48 (48) 




SpARCS J105111+581803 


1.035 


10:51:11.232 


58:18:03.128 


4 


194 


145 


34 (32) 




SpARCS J161641+554513 


1.156 


16:16:41.232 


55:45:25.708 


5 


242 


162 


46 (46) 


Demarco et al. (2010a) 


SpARCS J163435+402151 


1.177 


16:34:35.402 


40:21:51.588 


6 


205 


125 


50 (44) 


Muzzin et al. (2009) 


SpARCS J163852+403843 


1.196 


16:38:51.625 


40:38:42.893 


6 


189 


112 


44 (39) 


Muzzin et al. (2009) 


SpARCS J003550-431224[2,3l 


1.335 


00:35:49.700 


-43:12:24.160 


3 


91 


49 


20 


Wilson et al. (2009) 



^ The number of confirmed cluster members with LogM.^ /Mq > 9.3 is sliown in braclccts. 

This cluster is omitted in the current analysis because the usable wavelength coverage of the spectra is limited due to the high redsliift. 
Additional masks for this cluster have been observed and will be presented in the GCLASS catalog paper. 




Fig. 1. — gz[3.6/.tm] color images of the nine GCLASS clusters used in this analysis. The field of view of each image is 7' X 7', which is 
approximately the area covered by the spectroscopic observations and corresponds to 3.5 Mpc on a side at z = 1. Rich clusters at z ~ 1 
show a range of morphologies, from roughly spherical and centrally-concentrated with a clear central galaxy (e.g., SpARCS J003645-441050, 
z = 0.869), to asymmetric with filamentary-like structure and no clear central galaxy (e.g., SpARCS 1047374-574137, z = 0.956). 



6 



galaxies out to a radius of 3', which is ~ 1.5 Mpc at z ~ 
1. 

For all observations we used the OG610 filter which 
blocks light blueward of ^ 6IOOA observed frame. This 
made the spectra shorter so that two tiers of slits could 
be placed per mask. Each mask was observed with six ex- 
posures, each with 15 N&S cycles with 60s exposure time 
at each of the two nod positions. We nodded along the 
slit, which resulted in 30 min on-sky, per exposure, for a 
total of 3 hrs on-sky per mask. We used three positions 
for the central wavelength of the grating, 7380A, 7500A, 
and 7620A which dithers the spectra in the dispersion di- 
rection and effectively fills in the GMOS detector gaps. 
We also used the DTX offsetting in the spatial direction, 
which provides additional dithering, mitigating some of 
the noise from the N&S charge traps. 

Data reduction was performed using the IRAF GMOS 
package provided by the Gemini observatory. Each frame 
was subtracted using N&S darks. These frames follow 
the same shuffling pattern as the science data, but are 
taken with the shutter closed. Subtracting them removes 
the bias pattern, dark current, as well as charge that 
builds up and becomes trapped on certain pixels that do 
not shuffle the charge efficiently. Unlike typical spectro- 
scopic observations, we did not flat field the exposures 
before final combination because the GMOS chips, par- 
ticularly GMOS-S, suffer from significant fringing. For 
GMOS-S the fringing is ^ 10% peak-to-peak starting at 
A > 8OOOA but gets to as high as 55% peak-to-peak at 
A > 9OOOA. The N&S observation mode does a supe- 
rior job of removing the fringing in the spectra because 
the sky subtraction is determined using the same pixels 
that the object spectra is on (due to the on/off nodding 
pattern); however, the flat fields taken with an inter- 
nal lamp also have significant fringing. Attempting to 
correct for pixel-to-pixel sensitivity using the lamp ex- 
posures re-introduces pixel-to-pixel variations from the 
fringing that was removed with the N&S sky subtraction. 
The fringing is stronger than the pixel-to-pixel variations 
which are further mitigated by the dithering in both the 
spatial and dispersion directions, therefore the quality of 
data that is not flat fielded is substantially higher than 
that which is flat fielded. 

The six individual frames per mask were coadded 
into a final image of 2d spectra, and then Id spectra 
were extracted using the iGDDS softwar e developed for 
the G emini Deep Deep Survey (GDDS, Abraham et al.l 
l2004f l. Wavelength calibration is performed by identify- 
ing sky lines, and redshifts are obtained by comparison 
with templates available in the iGDDS software. Relative 
flux calibration was determined using long slit observa- 
tions of standard stars in each program. Comparison of 
the standard stars shows the relative flux calibration re- 
mained stable over the several years of observations. We 
apply the same flux calibration curve to all spectra, re- 
gardless of position, which assumes that the chromatic 
response of all detector pixels is the same. This may in- 
troduce errors in the continuum shape of the spectra if 
there is variance in the GMOS chromatic response that 
are position-dependent, although at present there is no 
evidence to suggest this is a significant issue. 

Based on the completed observations up to March 2011 
we created a vl.O spectroscopic catalog for the clusters. 
The vl.O catalog contains observations for 45/46 total 



masks in GCLASS. Table 1 shows the number of slits 
and successful redshifts in the field of each cluster in the 
sample, based on the vl.O catalogs. 

4. DETERMINATION OF GALAXY PROPERTIES 

4.1. Cluster Membership 

Galaxies are considered to be cluster members if their 
velocity relative to the cluster velocity in the rest-frame is 
Av < 1500 km s^^. We experimented with more sophis- 
ticated methods of defining membership based on frac- 
tional values of the cluster velocity dispersion, but found 
this made little difference to the final sample of clus- 
ter members. This is primarily because the clusters are 
well-defined in velocity space (see G. Wilson, in prepa- 
ration) with the vast majority of objects within 2a of 
line-of-sight velocity dispersion, values that are typically 
less than Av < 1500 km s~^. Therefore, in the interest 
of simplicity, we adopt this method for defining mem- 
bership. This definition is also conveniently similar to 
the line-of-sight cuts used for defining environment in 
field spectr oscopic surveys (e.g., 'KauffmanneiaU |2004; 
ICooper eFa l. 2005; Peng et al. 2010). Based on this cri- 
teria, the total number of cluster members across the 
sample of 9 clusters is 437. The number of members per 
cluster is listed in Table 1. 

4.2. Field Sample 

Throughout this paper we study the properties of 
galaxies as a function of environment, where environment 
is defined by the clustercentric radius of the galaxy. We 
have chosen to use physical clustercentric radius in kpc 
as our measure of the cluster environment, rather than 
rescale to the virial radius of each cluster. The clus- 
ters have a fairly modest range of velocity dispersions (a 
factor of ^ 3, G. Wilson et al., in preparation), so phys- 
ical clustercentric radii should scale roughly the same as 
virial radius for each cluster. A physical clustercentric 
radius is also appealing because it does not require the 
assumption of virial equilibrium, something which may 
not yet be achieved in young galaxy clusters. 

Clustercentric radius is a robust measure of environ- 
ment and correlates well with other environmental es- 
timators such as local density or Nth nearest neighbor; 
however, even at large clustercentric radius clusters con- 
tain a population of backsplash galaxies which have al- 
ready completed an orbit that brought them into or near 
the cluster core. Therefore, it is useful to define a con- 
trol population of "field" galaxies that have no associ- 
ation whatsoever with a rich galaxy cluster. The field 
population will be drawn from a wide range of environ- 
ments: from small groups to voids; however, as clusters 
grow they accrete galaxies from all environments, hence 
the field sample should be representative of the popula- 
tion of galaxies that are infalling into the clusters at that 
epoch. 

The GCLASS cluster sample spans 0.85 < z < 1.20 
so data from the survey itself can be used to define a 
population of field galaxies that have a similar redshift 
range as the clusters, as well as similar color/magnitude 
selection criteria. The very broad color prioritization in 
the mask design means that numerous field galaxies with 
similar colors as the cluster galaxies also have reliable 
spectroscopy. We consider galaxies at 0.85 < z < 1.20 
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Fig. 2. — Left Panel: Observed z' - 3.6/xm color vs. 3.6/im magnitude of galaxies in the field of the cluster SpARCS J105111+581803 
(rod points). Only galaxies with R < 550 kpc are plotted. The dashed lines denote the 5(7 completeness limits on the photometry. Galaxies 
identified as spectroscopic members are plotted with blue diamonds, and those that are confirmed foreground/background galaxies are 
plotted with grey diamonds. The shaded boxes represent the priority levels for targets in the mask design process. Right Panel: Same as 
left panel but for galaxies at R, < 1200 kpc. Overall, the spectroscopic completeness at the various priority levels is high, and the priority 
levels are well-tuned to obtaining the maximum number of cluster redshifts. 



that have Av > 2000 km s^^ from the chister velocity 
to not be associated with the cluster and are assigned 
to the field sample. This selection results in a sample of 
273 field galaxies with identical color/magnitude selec- 
tion criteria as the cluster galaxies. 

Given that massive clusters often form in overdense re- 
gions such as super clusters, there is some concern that 
the field sample will not be completely representative of a 
blank-field field sample, because line-of-sight groups may 
be more prevalent when looking towards a cluster, partic- 
ularly clusters selected as overdensities in an optical/IR 
survey. As a test, we drew an additional field sample of 
gala xies in the same redsh ift range from the GDDS sur- 
vey () Abraham et al.ll2004 l. The GDDS spectroscopy is 
taken with the same instrument as GCLASS and has the 
same resolution and uses the N&S mode, hence the spec- 
tra can be compared directly with the GCLASS spectra. 

The GDDS is a K-selected sample (K < 20.6) in four 
independent GMOS fields and has substantially deeper 
spectra than GCLASS (20-30 hrs on sky compared to 3 
hrs). In the 0.85 < z < 1.20 redshift range the observed 
K - 3.6/im colors of galaxies vary weakly with spectral 
type, hence it is straightforward to consider the GDDS 
sample in ter ms of a 3.6/xm-selected sa mple. Based on 
an average of (iBruzual fc Charlotl[2?)03l hereafter BC03) 
models with different star-formation histories and ages, 
we infer that at z = 1, the typical K-3.6/im color is K- 
3.6/im =1.2 and convert the K-band fiuxes of the GDDS 
galaxies to 3.6/xm fiuxes using this transformation. 

Once the K-band flux is converted to a 3.6/xm flux, 
the GDDS selection is 3.6^m < 19.4 (Vega magnitudes), 
which is nearly identical to the depth of the SWIRE ob- 
servations, 3.6/xm < 19.2 (Vega magnitudes). We com- 
pared the GCLASS field sample to the GDDS field sam- 
ple by performing a stacking analysis of the spectra (see 
§ 5). The equivalent width (EW) of the [Oil] emission 
line, and D„(4000) measures for the stacks are the same 
to within 5%, suggesting that the GCLASS field sample 
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• Constant SFR 
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Fig. 3. — Bottom Panel: Age tracks of models that show the 
M/L ratio in the observed 3.6/im band at z ~ 1 as a function of 
D„(4000). The age in Gyr (1 - 13) starts at the left side and is 
denoted by the circles on each track. Ages > 6 Gyr are forbidden by 
the age of the universe at z ~ 1 and are shown as open circles. Top 
Panel: The distribution of D„ (4000) values in the GCLASS cluster 
galaxies. The models show that galaxies can only obtain D„ (4000) 
> 1.7 by z ~ 1 if they have short star formation timescales and ages 
nearly as old as the universe. Even for the extreme star-formation 
histories considered here the range in M/Ls.g^m values at fixed 
D,i(4000) is at most a factor of a few. Thefore, the M/Ls.e^m for 
all galaxies is determined based on the t = 0.3 Gyr model. 

is representative of a blank-field field sample. Given the 
superior signal-to-noise ratio (S/N) of the GDDS spec- 
tra, we include the 87 GDDS galaxies at 0.85 < z < 1.20 
in our field sample, bringing the total number of galax- 
ies in the field sample to 360 galaxies. Colors, SFRs, and 
stellar masses for the GDDS galaxies are computed using 
the same methods as for the GCLASS galaxies. 

4.3. Stellar Masses 

Stellar masses are calculated for all galaxies by using 
the spectra to determine the stellar mass-to-light ratio 
(M(7i/L(n) in a manner similar to that developed by 
iKauffmann et al.l (I2003D for th e SDSS spectroscopic 
sample. iKauffmann et all (|2003f) used the BC03 models 
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to make synthetic spectra of galaxies with a range of 
star-formation histories. From these synthetic spectra 
they constructed a grid of observed z'-band stellar mass- 
to-light ratios {Mq/Lz') as a function of D(4000) and 
the H(5 absorption measure. Our method is similar, but 
with a few modifications given the different photometric 
bands in the surveys, and the fact that the typical S/N 
for GCLASS spectra is lower than that in the SDSS. 

We also use the BC03 stellar populatio n models with 
solar metallicity and a IChabrieii ([2003') initial mass 
function (IMF); however, we determine the observed 
3.6^m stellar mass-to-light ratios (MQ/La.g/im)- At 2; ~ 
1 this corresponds to the rest-frame H-band, which is 
similar to the z'-band for SDSS galaxies in that it has a 
weak dependence on galaxy spectral energy distribution 
( SEP) type a nd su ffers little from dust extinction. 

iMarastonI (|2005[ ) has released an alternative set of 
stellar population models which use a different treat- 
ment of the thermally-pulsating asymptotic giant branch 
stars (TP-AGB) than the BC03 models. Those models 
predict Mq/Lq in the rest- frame near- infrared that are 
approximately a factor of ~ 1.6 smaller than those from 
the BC03 models. At present it is unclear which of the 
available models provi des a better description of high- 
redshif t galaxies (e . g.. IWuvts et al.1 120071 : iMuzzin et al.l 
[2009at iKriek et al.l 120100 . so we have chosen to use 
the BC03 models to provide continuity with previous 
studies. We note that the choice of models does not 
affect our analysis in any significant way, because the 
difference in stellar masses between the BC03 and the 
Maraston models is a systematic offset of a factor of 
~ 1. 6, with little dependence on the stellar population 
(e.g., IWuvts et all 120071 : iMuzzin et aHl2009a[ ). 

The main differenc e between our method fo r determin- 
ing stellar mass and iKauffmann et al.l ()2003[ ) is that we 
use the D(4000) measures to estimate the Mo/La.e^m, 
but do not use the US measures. US is a relatively weak 
absorption feature with ~ l-SA EW for most galaxies. 
The S/N for the brighter galaxies in GCLASS is high 
and measuring H5 is trivial; however, this is not the 
case for the fainter spectra. The difficulty is that at the 
redshift range of the galaxy sample the rest-wavelength 
of US (4IO2A) corresponds to observed wavelengths of 
7590A to 9025A, a region that has numerous bright 
sky lines that at the low-resolution of the spectra cause 
significant noise. As we will show below, is an 
indicator of a young stellar population, but the fact 
that the universe is much younger at z = 1 than z = 
means that US is less important when determining 
Mq/L3 .6 /im- 

Like IKauffmann et all (|2003l ) we adopt the 
iBalogh et al.l (|1999D definition of D(4000) which is 
the ratio of in the flux density at 3850A-3950A to the 
flux density at 4000A-4100A, and hereafter refer to 
this as D„(4000). We parameterize the star- formation 
history of the galaxies using the standard declining 
exponential r-models of the form SFR cx e^*/'^, where t 
is the time since the onset of star formation, and r sets 
the star formation timescale. As an example. Figure 
3 shows the evolution of the observed Mo/La.g^m at 
z = 1 as a function of D„(4000) for models with a range 
of T. Each point on the tracks indicates an additional 
Gyr in age starting from the left side of the plot. Ages 



older than the age of the universe at z = 1 (t = 6 Gyr) 
are forbidden and are indicated by open circles. The 
measured D„(4000) values from the GCLASS spectra 
are shown as the histogram in the top panel of Figure 3. 

The arrow in Figure 3 shows the effect of 1 magnitude 
of V-band extincti on (A„ = 1.0) on the D„ (4000) of the 
models assuming a ICalzetti et all (|2000[ ) dust law. This 
weak dependence on dust illustrates an advantage of 
using D„(4000) as an indicator of the M0/L0 compared 
to SED-fitting. This weak dependence occurs because of 
the narrow wavelength range used to measure D„(4000), 
compared to the much broader dust extinction curve. Of 
course, the disadvatange of using D„(4000) as compared 
to SED-fitting is that it does not provide additional 
information on the star-formation history of the galaxy, 
i.e., which r- model should be used. 

Figure 3 shows the maximum range in MQ/Ls.g^m 
between the most extreme models (the youngest con- 
stant star formation vs. the oldest single-burst) is 
a factor of ^ 13. The age of the universe does not 
allow for populations with t > 6 Gyr, and therefore for 
galaxies in the GCLASS sample the maximum range 
in allowed MQ/Ls.g^m is only half as large, a factor of 
~ 7. This is still substantial; however, an additional 
restriction on the MQ/La.g/^m is the fact that not all 
star-formation histories are consistent with the data. 
Models with r > 1 Gyr do not provide sufficiently strong 
D„(4000)'s to match many of the cluster galaxies, even 
if the stellar populations are as old as the age of the 
universe. Many of the D„(4000)'s are so strong they can 
only be obtained with relatively short star formation 
timescales (r < 0.3 Gyr), and old ages. For the 94% of 
galaxies with D„(4000) > 1.2, the maximum difference 
in allowed Mo/Ls.g/im at fixed D„(4000) is only a 
factor of ^ 4. This factor is the most extreme range 
in allowed Mo/La.g^m across all models and suggests 
that the stellar masses would be reasonably accurate 
even without any additional information on the star 
formation history from the D„(4000). 

Simi lar to other high-redshift st el lar mass studies 
(e.g., iForster Schreiber et al.l 120041 : iMarchesini et al.l 
l2009f ) we adopt the r = 0.3 Gyr model as our default 
model and thereafter determine the Mo/La.g/im based 
on the D„(4000). This approach may underestimate 
the M0/L3.6/j,m for populations with D„(4000) < 1.2 
by as much as a factor of ~ 2, but should work well 
for the majority of the galaxies in the sample. Given 
the other systemati c uncertainties involved in stellar 
mass modeling (e.g.,lMaraston et al.ll2006l : IMuzzin et al.l 



I2009at [CM7rov et al.l 



2OO90 . and that our primary use of 



the stellar masses is to rank-order galaxies in increasing 
mass, this should not affect our analysis. We note that 
it may cause systematic offsets when comparing our 
stellar masses of galaxies with weak D„(4000) to those 
determined by other studies with the more typical SED 
fitting method, and advise caution if doing so. 

The full spectroscopic sample of galaxies consists of 
437 cluster members and 360 field galaxies. Of the 
437 cluster members, 9 of the galaxies are not detected 
at 3.6/^m, and of the 360 field galaxies, 26 are not 
detected at 3. 6/im. Based on the SWIRE 3. 6/im flux 
limit we estimate that these galaxies have LogM,/MQ 
< 9.3. For these galaxies we do not attempt to calculate 
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stellar masses using the D„(4000) models and their 
z'-band magnitudes. At z ^ 1 the z'-band corresponds 
to rest-frame B-band. The Mq/L^/ at z = 1 have a 
substantially larger range at fixed D„(4000) than the 
M0/L3.6/im, and also a strong dependence on dust 
content and hence are uncertain. For the remainder 
of this analysis we remove these galaxies from the 
sample and only consider galaxies above a limiting 
mass of LogM,/M0 > 9.3. This requires removing and 
additional 16 cluster galaxies and 40 field galaxies that 
have 3.6/im detections, but LogM,/M0 < 9.3; however, 
it provides a firm 3.6/im-limited sample, which can 
be corrected for completeness to a stellar-mass-limited 
sample (see § 4.6). Once the limit of LogM,/M0 > 9.3 
has been applied, the final sample of cluster members 
contains 412 galaxies and the final field sample contains 
294 galaxies. The number of cluster members in each 
cluster above the mass limit is indicated in Table 1. 

4.4. Star Formation Rates 
Star formation rates are determined using the lumi- 
nosity of the [Oil] emission line (L3727A line)' '^hich 
has been determined from the aperture-corrected line 
flux (F3727A,iine)- The F3727A,iine is itself derived from 
the measurement of the EW of the [Oil] line. The 
first step in determining the EW([OII]) is to fit the 
emission line to a Gaussian in the rest-frame. The 
continuum level ^^^27}^ cont then determined using 
the region of the spectrum within 50A on both the 
red and blue side of the emission line. The EW([OII]) 
is then converted to an aperture-corrected line fiux 
using the total 3.6/im flux. We determine the ra- 
tio between the continuum at rest-frame 3727A and 
observed-frame 3.6/im ^^•j27k cont 

/F3 .6/im,obs)model 

from the BC03 model implied by the D„(4000) mea- 
surement (see § 4.3). This is measured for each 
galaxy individually because the rest-frame probed by 
3.6/im is slightly different depending on the redshift 
of the galaxy. The [Oil] line flux is then ^^■j2-jk line 

— F3,6/im X (F3727A cont 

/F3.6Mm,obs)modol XEW([0II]), 

where F3 is the observed total 3.6/im flux of the 
galaxy 

Once the ^^^271^ line is determined we convert this to 
a L3727A line ^sing the spectroscopic redshift and apply 
the mass-dependent [OII] -SFR relation determ ined by 
iGilbank et all (|201^. The' Gilbank et al.l (|2010l ) relation 
was empirically calibrated using Balmer-decrement 
corrected Ha SFRs in the SDSS. This relation therefore 
encodes the extinction correction and metallicity de- 
pendence of the [Oil] line implicitly by assuming these 
quantities are correlated with stellar mass. As shown 
in Figure 3 of Gilbank et al., the mass-dependence of 
the conversion is weak for galaxies with LogM*/M0 
> 10.5, which comprise much of our sample, but gets 
increasingly larger for galaxies with LogM*/M0 < 10.5. 
We note that although the mass-dependence of the 
correction is weak at LogM»/M0 > 10.5, the correction 
itself is significant, approximately a factor of 3 larger 
than a constant-luminosity [Oil] SFR. 

4.5. Rest-frame U-Bghort Colors 



Synthetic rest-frame colors for the galaxies are cal- 
culated by observing the flux-calibrated spectra in the 
rest-frame with the appropriate filter response functions. 
The wavelength coverage of the spectra in the observed 
frame is 6IOOA - 9500A, which corresponds to rest-frame 
33OOA - 5140A for the lowest-redshift galaxies (z = 0.85), 
and 27OOA - 4320A for the highest-redshift galaxies (z = 
1.20). The wavelength coverage is sufficient to compute 
a synthetic U-band flux for galaxies at all redshifts; how- 
ever, the B-band is only completely within the spectral 
window for the lowest-redshift galaxies. The peak trans- 
mission for the B-band is ^ 4100A; however, it has a long 
tail of transmission redward that does not decrease to < 
10% of maximum transmission until ^ 5250A. Therefore, 
in order to compute a rest-frame color, we define a modi- 
fied B-band filter where we set the transmission redward 
of 4450 A equal to zero. This synthetic filter, which here- 
after we refer to as "B^^ort" allows us to compute a U 
- B color for all but the highest redshift galaxies in the 
sample, but still retains the general characteristics of the 
standard B-band. 

We tested the effect of shortening the B-band on galaxy 
colors by comparing the U - B colors of synthetic spec- 
tra from the BC03 models with a range of colors to the 
U - ^short colors for the same models. The transfor- 
mation changes as a function of galaxy color, but for 
galaxies with U - ^short > 0.6, the majority of objects in 
the GCLASS sample, the transformation is U - B ~ U - 
Bs/tort + 0.15. 

4.6. Spectroscopic Completeness 

The completeness of a spectroscopic sample depends 
both on the sampling frequency of the targets, as well as 
redshift success rate. In turn, both of these depend on 
observable quantities of the galaxies such as flux, color, 
position, redshift, and the presence/lack of observable 
emission lines. The sampling frequency depends on ob- 
servables because galaxies are prioritized in the mask de- 
sign process based on their properties. The redshift suc- 
cess rate depends on observables, particularly flux and 
redshift, because these quantities increase or decrease 
the confidence of the redshift determination. Brighter 
targets provide higher S/N spectra, and galaxies at par- 
ticular redshifts have their prominent spectral features 
within the observed bandpass. If the primary source of 
incompleteness is the sampling frequency, not the red- 
shift success rate, it is straightforward to compute spec- 
troscopic completeness corrections. This condition can 
usually be met by limiting the overall sample to within 
a fiux and redshift limit where the redshift success rate 
is high. 

Due to the relatively deep exposures and good N&S 
sky subtraction the redshift success rate for all slits, re- 
gardless of brightness and redshift is quite good (71%). 
In particular, the success rate for the primary priorities, 
PI - PIG, ranged between a minimum of 70% (P7) to a 
maximum of 96% (PI). Even the bright filler slits (Pll) 
had a good redshift success rate of 69%. The major- 
ity of slits where no redshift could be obtained were the 
faint filler slits (P12) which had a redshift success rate 
of only 21%. The majority of the P12 objects fall below 
the fiux/stellar mass limit that we have adopted for the 
sample (LogMH./M0 > 9.3, see below), and only 1 was a 
cluster member. We examined the spectra of PI - PIO 
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Fig. 4. — Left Panel: Spectroscopic completeness fraction of cluster galaxies as a function of clustercentric radius and stellar mass. The 
completeness is calculated based on the color and magnitude space occupied by cluster members in each cluster, and the relative fraction 
of galaxies in the color space that have redshifts (see text). Spectroscopic completeness is highest for massive galaxies and galaxies in the 
cluster core. Right panel: Spectroscopic completeness as a function of radius and stellar mass normalized to the completeness for massive 
galaxies. This shows that the primary completeness bias is stellar mass; however, this bias is similar at all radii. 



galaxies in order to ascertain why the redshift success 
rate was less than 100%. For many PI - PIO galaxies 
where no redshift was determined the spectrum did have 
clearly detected continuum emission; however, there were 
no identifiable spectral features. It is likely that most of 
these targets are either galaxies at z < 0.5 where the 
band-limiting filter at 6IOOA makes it difficult to iden- 
tify features, or galaxies at z > 1.5 where there are few 
strong spectral features in the observed bandpass. There- 
fore, provided we limit the sample to the redshift range 
where the spectral features are in the observed bandpass 
(0.85 < z < 1.20), and to galaxies above the flux/stellar 
mass limit required for a successful redshift (LogM^/M© 
> 9.3), the completeness corrections can be calculated 
under the assumption that the sampling frequency is the 
dominant source of incompleteness. 

Despite the fact that the completeness depends on ob- 
servable quantities, the spectroscopic completeness cor- 
rections need to be calculated as a function of the physi- 
cal quantities that are used to define the galaxy sample. 
For the current study, where we examine the properties of 
galaxies as a function of stellar mass and environment, 
we calculate completeness correction in terms of those 
two quantities. 

As was shown in § 4.3, the stellar mass of galaxies 
is closely related to their 3.6^m flux, the primary se- 
lection criteria for spectroscopic targets. In principle, 
this means that defining the completeness above a stel- 
lar mass limit is nearly analogous to defining it above a 
3.6/im flux limit. 

For each cluster, the quantity of interest is the spec- 
troscopic completeness for cluster galaxies, not all galax- 
ies in the field of view. This means we are not in- 
terested in the completeness of a flux-limited sample, 
but a flux-limited sample within a particular region of 
color-magnitude space, i.e., the color- magnitude space 
occupied by z ^ 1 cluster members. This region of 
color-magnitude space will be different from cluster-to- 
cluster because they are at a range of redshifts and the k- 
correction affects the color of the observed red-sequence 
and blue cloud. 

The location of the cluster in color-magnitude space 
is determined by examining the location of all con- 
flrmed cluster members in the z' - 3.6/im vs. 3.6/im 



color-magnitude diagram. We then define the color- 
magnitude space inhabited by cluster galaxies using sev- 
eral adaptive-size boxes that are adjusted to cover the 
color space occupied by all cluster members. The com- 
pleteness for cluster members as a function of 3.6/fm flux 
within that region of color space is then defined as 



fclus 



Ngal{Ncluster,z / (N field, z + Nduster.z)) ' 



(1) 



where 'Lciuste.r,z is the fraction of cluster members in that 
flux range that have a spectroscopic redshift, ^ciuster,z 
is the total number of cluster members in the color- 
magnitude space with a redshift, ^fieid,z is the total 
number of field galaxies in the color-magnitude space 
with a spectroscopic redshift, and Ng^i is the total num- 
ber of galaxies with or without a redshift in the color- 
magnitude space. 

Once the completeness in this color space is computed, 
we convert this to a completeness as a function of stellar 
mass using a fit to the correlation between 3.6/im and 
stellar mass. The correlation between these two param- 
eters is linear in log space and has an rms scatter of 0.15 
dex. For the final completeness calculations we combine 
the data from all 9 clusters into an ensemble cluster and 
consider the completeness in several bins of stellar mass, 
and clustercentric radius. The completeness curves for 
galaxies of different stellar masses are plotted as a func- 
tion of radius in the left panels of Figure 4. 

Figure 4 shows clearly that the spectroscopic complete- 
ness is a function of both stellar mass, and clustercentric 
radius. The completeness is highest for massive galaxies 
in the cluster cores, and lowest for lower-mass galaxies 
in the cluster outskirts. Both of these incompleteness ef- 
fects are straightforward to understand. 

The increasing completeness with increasing stellar 
mass is primarily a product of the 3.6/im flux prioriti- 
zation and the fact that bright galaxies have the highest 
S/N spectra, hence obtaining redshifts is easiest. The de- 
creasing completeness with increasing clustercentric ra- 
dius is primarily an area effect. Due to the microslits 
available with the band shuffle spectroscopy mode, the 
total number of redshifts as a function of clustercentric 
radius is roughly constant; however, the area grows as 
7rr^, and hence obtaining redshifts for the same fraction 
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of galaxies becomes increasing difficult at larger radii. 

In the right panel of Figure 4 we plot the spectro- 
scopic completeness of galaxies as function of cluster- 
centric radius, scaled to the same completeness as a func- 
tion of clustercentric radius for the most massive galax- 
ies. These curves are not perfectly flat but have only 
mild slopes, demonstrating that the main incompleteness 
in GCLASS is the stellar mass incompleteness, and that 
the stellar mass incompleteness is nearly independent of 
environment. 

The main focus of this paper is to examine the prop- 
erties of galaxies at fixed stellar mass, and fixed environ- 
ment (see § 6). As the right panel of Figure 4 shows, 
the sampling for galaxies of a given stellar mass rela- 
tive to other stellar masses is similar in all environments. 
Therefore, completeness corrections are not required in 
that analysis. In § 5, where we examine the properties 
of galaxies for a mass-limited sample, we do employ the 
completeness corrections as a function of stellar mass and 
environment to correct for those biases. 

5. RESULTS 

Most of the analysis presented in this paper is based on 
measurements made from stacking the spectra. The main 
advantage of stacking is that it increases the S/N which 
improves the precision of the measurement of the mean 
properties of galaxies. Of course, the tradeoff in bet- 
ter measurements of the mean is that information about 
the distribution in properties is lost. Given that there 
is a large range in S/N for the spectra, the distribution 
of properties is not always easy to interpret, and there- 
fore we have concentrated on robust measurements of the 
mean properties of the galaxies. 

When stacking, each spectrum is normalized using the 
median flux in the range 4050A-4100A, rest-frame. This 
part of the spectrum is redward of the 4000A break, but 
does not contain substantial absorption or emission lines. 
After normalization, each galaxy is weighted by its po- 
sition and stellar mass completeness weight (see § 4.6). 
The rest-frame wavelength coverage of the spectra de- 
pends on the redshift of the galaxy, and therefore addi- 
tional weighting is performed based on how many galax- 
ies contribute to the stack on the red and blue ends. Once 
all weighting has been performed, the individual spectra 
are combined into a mean stacked spectrum. 

Throughout the analysis, the uncertainties in proper- 
ties measured from the spectral stacks (e.g., SFR and 
D„(4000)) are determined using 300 bootstrap resam- 
plings of the data. Many of the stacks contain > 40 
galaxies and hence the random errors in most properties 
are small, typically < 1%. In these cases systematic er- 
rors are almost certainly the dominant source of error. 
The systematics of greatest concern are the models em- 
ployed in the determination of physical parameters, such 
as the conversion of [Oil] line fluxes to SFRs, or the de- 
termination of stellar masses from synthetic BC03 spec- 
tra. At present, these types of systematics are difficult to 
quantify, and indeed are probably the dominant source of 
error in most studies that determine physical properties 
of galaxies from observables. Therefore, throughout this 
paper the quoted errors are always the bootstrap random 
errors, but we note that systematic errors are likely to 
be substantially larger and advise caution if comparing 
quantities such as SFRs and stellar masses determined 



in this paper to quantities determined in other papers 
using difii'erent model assumptions. 

5.1. Star Formation Rates as a Function of 
Environment and Stellar Mass 

Here we investigate the relationship between the SFRs 
and SSFRs of galaxies with their environment and stel- 
lar mass. In the left panel of Figure; 5 we plot the 
completeness-corrected mean stacked spectra of galax- 
ies with LogM*/M0 > 9.3 as a function of increasing 
clustercentric radius culminating with the field sample. 

Without making any quantitative measurements it is 
immediately clear from the left panel of Figiu-e 5 that the 
average stellar populations of galaxies are dramatically 
different depending on whether they are located in the 
field, the cluster outskirts, or the cluster core. Remark- 
ably, it appears that the transition in stellar popiilations 
as a function of clustercentric radius is quite smooth and 
that there is no preferred location where the properties 
of galaxies change distinctly. 

Galaxies in the cluster cores (R < 200 kpc) have clear 
4000A breaks and show the Mgl (3830A) and G-band 
(4304A) absorption features. Cluster core galaxies also 
exhibit weak [Oil] emission. All of these features are 
indications of an evolved stellar population with little 
ongoing star formation. 

The galaxies in the field sample are considerably dif- 
ferent. They show the Balmer series absorption lines and 
have a Balmer-break rather than a 4000A break. They 
also show evidence of [Nelll] (3869A) emission and have 
strong [Oil] emission, all of which are indications of a 
relatively young stellar population with ongoing star for- 
mation and possible AGN activity. 

The SSFRs and SFRs of the stacks are determined 
from the measured EW([OII]) using the method de- 
scribed in § 4.4 and are plotted in the left panels of 
Figure 6, as well as listed in Tables 2 and 3. There is 
a 1.33 ± 0.09 dex decline in the SSFR of galaxies as 
their environment changes from the field into the cores 
of rich clusters, making it clear that the average star 
formation properties of galaxies are strongly correlated 
with their environment at z ~ 1. Figure 6 shows that 
not only does the mean SSFR decrease with increasing 
galaxy density, but that the mean SFR of galaxies also 
declines from the field toward the cluster core. This con- 
firms that both the SSFR-density relation and the SFR- 
density relation are already in place at z ~ 1, at least 
in the highest-density environments. Interestingly, the 
SFR-density relation only declines by 0.67 ± 0.03 dex 
between the cluster core and field, not nearly as much as 
the SSFR-density relation. This suggests that the mean 
stellar masses of cluster and field galaxies are not the 
same, and may drive some of these correlations, an issue 
we will return to in §5.4. 

The trends in the SSFRs and SFRs with environment 
are quite clear; however, we know that these properties 
also correlate with galaxy stellar mass. In the right pan- 
els of Figure 5 we plot the mean stacked spectra of the 
cluster galaxies in several bins of stellar mass. Figure 
5 shows that the same trends that exist between stellar 
populations and increasing environmental density also 
exist with increasing stellar mass. More massive galax- 
ies appear to be a more evolved population with little 
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Fig. 5. — Left Panel: Completeness-corrected mean stacked spectra of cluster galaxies with LogM»/M0 > 9.3 as a function of clustercentric 
radius. Field galaxies are taken from both GCLASS and the GDDS (see text). There is a clear change in the spectral properties of galaxies 
with clustercentric radius, moving from 4000A break dominated systems with a weak [Oil] 3727A in the core to Balmer-break dominated 
systems with a strong [Oil] 3727A in the field. Right Panel; Completeness-corrected mean stacked spectra of cluster galaxies (R < 1.5 
Mpc) as a function of stellar mass. The change in properties with increasing stellar mass appears to mirror the change in properties with 
increasing environmental density. 



TABLE 2 

Galaxy Parameters for the Mass-limited Sample as a Function of Environment 
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TABLE 3 

Galaxy Parameters for the Mass-limited Sample as a Function of Stellar Mass 


Log(Mass) 
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ongoing star formation, whereas lower-mass galaxies are 
younger with clear signs of active star formation. 

In the right panels of Figure 6 we plot the inferred 
SSFR and SFR from these stacks as a function of stellar 
mass. For reference we also performed a stacking analy- 
sis of the field sample as a function of stellar mass and 
also show those SSFRs and SFRs in Figure 6. As with 
environment, there is a clear trend of decreasing SSFR 
with increasing galaxy stellar mass both for galaxies in 
clusters as well as in the field. Interestingly, the middle 
panel of Figure 6 shows that the mean SFR of the galax- 
ies is a fairly weak function of stellar mass. This indepen- 
dence of SFR on galaxy mass is also seen in the HiZEL S 
narrow-band Ha survey at z ~ 0.84 f Sobral et aLll201lD . 
and is most likely a co-incidence of the overall process of 
galaxy quenching as a function of stellar mass and red- 
shift and the particular redshift range considered here. 

The results shown in Figure 6 are in good agreement 
with most of the previous work on z ^ 1 clusters, but do 
not agree as well with the z ~ 1 field galaxy studies. The 
1.33 ± 0.09 dex decline in SSFR seen with e nviro nment 
is similar to the study of lPatel et all ()2009l |2011D . who 
measured the SSFRs and SFRs of galaxies around a mas- 
sive galaxy cl uster at z = 0.83 u sing low-resolution grism 
spectroscopy. iPatel et all ()2011|) measured a 0.85 dex de- 
cline in the SSFR for all galaxies in their mass-limited 
sample (LogM^/M© > 10.25) from the low-density field 
to the core of their rich cluster. 

The trend of decreasing SSFR with increasing galaxy 
density is qualitatively similar to the DEEP2 study of 
ICooper eF al. (2 00^; h o wever , quantitatively it is quite 
different. Co oper et al.l (|2008[ ) measured only a 0.14 dex 
decline in SSFR between the highest and lowest density 
environments found in DEEP2 at z ~ 1 , a negligible envi- 
ronmental effect. As a comparison we plot the mean field 
value from the DEEP2 study beside our field sample in 
Figure 6. Considering the differences between GCLASS 
and DEEP2: 3.6um-spectroscopic selection vs. R-band, 
and a different [Oil] extinction correction, these values 
are in reasonably good agreement. This suggests that 
the discrepancy between our studies is not from different 
methods of computing SSFRs. It may be that the effects 
of environment on the galaxy population at z ~ 1 are 
relatively minor in all environments except for the rich- 
est galaxy clusters, so mething that has been suggested 
by other authors fe.g.. ISobral et al.l l2011h. 

5.2. The Fraction of Star-Forming and 
N on- Star- Forming Galaxies as a Function of 
Environment and Stellar Mass 

In the bottom panels of Figure 6 we plot the fraction 
of galaxies with detectable [Oil] emission as a function 
of environment and stellar mass, and hereafter refer to 
these galaxies as star-forming galaxies. The star-forming 
galaxies have been identified by examination of individ- 
ual object spectra. We have tested our ability to iden- 
tify [Oil] emission at different S/N levels by degrading 
the S/N of the best spectra with [Oil] emission. These 
tests show that [Oil] emission can be identified down to 
lA equivalent width for the highest S/N noise spec- 
tra, and down to ^ 3A equivalent width for the low- 
est S/N spectra. For the typical galaxy in our sample 
an EW([OII]) of 3A corresponds to a SSFR of - 5 x 



10^^^ yr^^. Formally this means that some galaxies 
we consider quiescent may not be completely quenched; 
however, given that the SSFRs of galaxies decline quite 
sharply at z < 1, galaxies with SSFRs this low are un- 
likely to add significantly more stellar mass through star 
formation in the future. We tested if we may have mis- 
classified a substantial population of star-forming galax- 
ies with low SSFRs as quiescent by stacking the spectra 
of the quiescent galaxies. We find no evidence for a de- 
tection of [Oil] in the quiescent stacks at a level of < 
lA, suggesting that > 90% of the galaxies we classify as 
quiescent have SSFRs < 5 x 10^^^ yr^"'^. 

Figure 6 shows there is a strong correlation between 
the fraction of star-forming galaxies with both environ- 
ment and stellar mass. The trends with environment 
are the strongest, with 82% ± 7% of field galaxies be- 
ing star forming, but only 22% ± 5% of galaxies in the 
cores of clusters showing star formation. The correlation 
between is f and environment ha s been observed ubiqui- 
tously in the local u niverse (e-g- . Kauffma nn et al.l[2004t 
iBalogh et aLll2004bt iPeng et aLlfeOlO ) and our data con- 
firms it is still in place at z ~ 1. O ur isF as a functio n 
of environment also agrees well with iPatel et al.l (|2011t ) , 
who found 79% ± 4% in low density enviro nments and 
32% ± 3% in the cores of rich clusters. The 'Pate l et al.l 
(j2011t ) galaxies are classified based on their location in 
the rest-frame U-V vs. V-J color-color diagram, not [Oil] 
emission, but the consistency between the measurements 
suggest both methods identify similar galaxies. 

There is also a change in the isF as a function of galaxy 
stellar mass both for field and cluster galaxies. For clus- 
ter galaxies, the percentage of galaxies in the highest stel- 
lar mass range (LogM^/M© > 11.0) that are star forming 
is 27% ± 7%, whereas in the lowest stellar mass range 
(9.3 < LogM,/M0 < 10.0), 56% ± 13% of the galax- 
ies are star forming. For field galaxies the corresponding 
numbers are 58% ± 18% and 96% ± 13%, respectively. 
This shows that a significant portion of the correlation 
between both SSFR and SFR with stellar mass and en- 
vironment must be driven by the changing fraction of 
star-forming galaxies with these parameters. This is an 
important distinction that is different from an overall de- 
cline in the SSFRs and SFRs of the star-forming galax- 
ies. In § 6 we will examine the mean properties of star- 
forming galaxies alone to ascertain what portion of the 
decline in SSFR and SFR with environment and stellar 
mass are caused by the change in fg p compared to a de- 
cline in the SSFRs and SFRs of the star-forming galaxies. 

5.3. Dn(4000) and U-Bshort colors as a function of 
Environment and Stellar Mass 

In the left panel of of Figure 7 we plot the measured 
D„(4000)s and \J-B short colors of the stacked spectra as 
a function of clustercentric radius. There is a correlation 
of D„(4000) with clustercentric radius, with D„(4000) 
changing from a value of 1.213 ± 0.003 in the field to 
1.606 ± 0.019 in the cluster cores. The D„(4000) is pro- 
portional to the luminosity- weighted age and metallicity 
of a galaxy, hence we can infer that the galaxies in the 
cluster cores have on average older and/or more metal- 
rich stellar populations than those in the field. This 
inference is fully consistent with the spectral features 
observed in the stacked spectra (Figure 5), where field 



14 




Fig. 6. — Left Panels: Log(SSFR), Log(SFR), and Igp for galaxies with LogM,/M0 > 9.3 as a function of clustercentric radius. The 
SSFR and SFR are determined from the stacked spectra in Figure 5. Error bars are calculated from 300 bootstrap resampling of the stacks 
but in same cases are smaller than the data points. The median value of the DEEP2 field sample which has as a similar redshift to the 
GCLASS sample is shown for comparison. At ^ ~ 1 there is clearly a strong SSFR-density and fgij-density relation in place. There is 
also a SFR-density relation, but it is weaker than the other two correlations. Right Panels: Same as left panels but plotted as a function 
of stellar mass both for cluster galaxies (magenta) and field galaxies (green). The cluster galaxies display the typical "downsizing" trends 
with more massive galaxies being less frequently star forming and having lower SSFRs. 




Fig. 7. — Left Panels: D„(4000) and rest-frame U - Tigf^^rt color for galaxies with LogMt/M© > 9.3 as a function of clustercentric radius. 
The Dn(4000) and U-Ba^o^^ color are determined from the stacked spectra in Figure 5. Error bars are calculated from 300 bootstrap 
resampling of the stacks but in same cases are smaller than the data points. There is a clear Dn (4000)-density and color-density relation in 
place at 2 ~ 1. Right Panels: Same as left panels but plotted as a function of stellar mass. Similar to the local universe both parameters 
are a strong function of galaxy stellar mass. 
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galaxies are dominated by B aimer lines, and cluster core 
galaxies have the metal lines commonly found in older 
stellar populations. 

There is also clear evidence for a color-density rela- 
tion, with galaxies in the cluster cores having \J-Bshort 
colors 0.38 ± 0.02 mag redder than those in the field. 
This is similar t o the color-density r elation seen at z ^ 
1 in DEEP2 bv ICooper etall (|20'06[) . They found that 
the typical U-B color of galaxies changes by ~ 0.8 mag 
between the highest and lowest density environments. It 
is surprising that the DEEP2 color-density relation is 
stronger than the color-density relation measured in the 
GCLASS clusters, but that the DEEP2 SSFR-density re- 
lation is substantially weaker than that measured in the 
GCLASS clusters (§ 5.1). Nonetheless, both studies do 
confirm a significant color-density relation at z ~ 1. 

In the right panel of Figure 7 we plot the measured 
D„(4000) and V-Bghort color for galaxies in different stel- 
lar mass bins. The trends in these parameters with stellar 
mass are remarkably similar to the trends with environ- 
ment, and it is clear that there is also a D„(4000)-stellar 
mass and color-stellar mass relation at z ~ 1 as well. 

5.4. The Dependence of Stellar Mass on Environment 

In the previous sections we have shown that the cor- 
relations between galaxy properties such as SSFR, SFR, 
D„(4000), V-Bshort, and isr with both their environ- 
ment and their stellar mass at z ~ 1 are quite strong. It 
also appears that the effect of increasing stellar mass on 
the properties of galaxies is identical to the effect from in- 
creasing environmental density. Indeed, the correlations 
of all parameters with stellar mass and environment are 
even quantitatively quite similar - although it is impor- 
tant to note that this may be a coincidence of the partic- 
ular range of environments and stellar masses considered 
in this study. 

Nonetheless, from the analysis up to this point we 
can safely conclude that the properties of galaxies are 
strongly correlated with both their stellar mass and en- 
vironment at z ~ 1. However, we cannot yet take the 
next step and conclude how causal those correlations 
are. If the stellar mass of galaxies and their environ- 
ment are themselves correlated it may be that only one 
of them is causally implicated in galaxy evolution and 
that the correlation between the other parameter and 
galaxy properties is simply a result of a more fundamen- 
tal stellar mass vs. environment correlation. Indeed, 
a correlation between galaxy stellar mass and environ- 
ment is observed by most stud i es of galaxies in th e lo- 
cal fe.g.. lKauffmann et al.l l2004':'Baldrv et al."2006'). and 
high-redshift universe (e.g., 'Bolzonella et al. 2010). 

In Figure 8 we plot the stellar mass of galaxies as a 
function of increasing clustercentric radius from the clus- 
ter cores to the field sample. Points in grey are galaxies 
with spectroscopic redshifts, but that fall below the mass 
limit of LogM,/Mo > 9.3, and are shown to illustrate the 
completeness. Galaxies detected in the z'-band but not 
at 3.6/xm are plotted as arrows. The vertical dashed line 
shows twice the FWHM of the IRAC PSF, and indicates 
where there may be incompleteness in lower-mass galax- 
ies due to blending with the central galaxy. The yellow 
line in Figure 8 shows the running mean determined by 
the 10 nearest galaxies, and does not include the grey 
points. We emphasize that the yellow line represents the 
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Fig. 8. — Stellar mass of galaxies as a function of clustercentric 
radius. Points in grey are galaxies that fall below the LogM*/M0 
> 9.3 completeness limit. The running mean in bins of 10 galaxies 
is shown as the yellow line. Field galaxies are shown at the right 
and have been offset randomly along the x-axis for clarity. In- 
cluding the central galaxies, there is evidence for mass segregation 
in the clusters with central galaxies being an order of magnitude 
more massive than satellites on the outskirts. Considering only the 
satellites, the mass segregation is weaker. On average, satellites in 
the cluster core are twice as massive as satellites on the cluster 
outskirts. In this plot the completeness is a strong function of 
both stellar mass and radius in that it is higher for massive galax- 
ies and galaxies in the core; however, the relative completeness in 
stellar mass is nearly independent of radius (see § 4.6) so the mass 
segregation is a real and not a selection effect. 



mean stellar mass of galaxies with LogM^/M© > 9.3, not 
the mean stellar mass of all galaxies in z ~ 1 clusters. 

Figure 8 shows that the mean stellar mass of galaxies 
decreases by ~ 1 dex from the cluster cores out to the 
field. This decline is substantial; however, it is driven pri- 
marily by the extremely massive cluster central galaxies 
which have few counterparts at larger radii. Ignoring 
the centrals and considering only satellites at R > 0.1 
Mpc shows that there is still a decrease in the mean stel- 
lar mass with environment, but it is much less, approx- 
imately a factor of 2. This covariance between mean 
stellar mass and environment demonstrates that in or- 
der to understand which parameter is most important in 
shaping the properties of galaxies, we must look at envi- 
ronmental effects at fixed stellar mass, as well as stellar 
mass effects at fixed environment. 

One concern with Figure 8 is that the stellar-mass in- 
completeness in the GCLASS spectroscopic sample may 
create a correlation between environment and stellar 
mass that does not exist. As we showed in § 4.6, lower- 
mass galaxies are under-represented by a factor of '~ 2-3 
compared to the most massive galaxies. Despite this, 
the large number of masks per cluster and ability to 
pack many slits in the cluster core with the band-shuffle 
mode means that for each cluster we are able to sample 
the same galaxies equally well in all parts of the clus- 
ter, and the stellar mass selection effects are basically 
independent of clustercentric radius (see Figure 4). In- 
deed, if anything, lower-mass galaxies are slightly under- 
represented at large radius. This means that not only 
is the trend of decreasing galaxy stellar mass with in- 
creasing clustercentric radius seen in Figure 8 is not a 
selection effect, but it may be even more significant than 
we have inferred. 
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TABLE 4 

Galaxy Parameters for Star, Forming Galaxies as a 

FrNCTiox OF Mass and Exmronmfnt 



Radius fgF Log(SSFR) D„(4000) Nstack 

Mpc yr~^ 

(1) (2) (3) (4) (5) 



Log(M,/M0) > 10.7 



0.1 
0.5 
1.0 
Field 


20+° °'^ 
0.3710;11 


-q 24+0-07 
-9 12+0-10 
-9 24+0-08 
-9 20+0-02 


1 624+°-°2* 
^•"^^-0.028 

1 542+°-°"^ 

1 624+0-028 
^•"^^-0.028 
1 519+0-010 


14 
10 
15 
43 


10.0 < Log(M*/M0) < 10.7 


0.1 
0.5 
1.0 

Field 


19+005 
26+0 04 
0.5410.10 

0.85l°-ii 


-8.90tVsl 

-8.r2l°°l 

o 07+0. 03 


1 256+0-0-1-1 
1 Qqo+0.012 
1-3'3^_0.012 
1 097+0.007 
1-^^' -0.007 

1 r,iQ + 0.003 


6 
32 

34 
99 




9.3 < Log(M*/M0 


) < 10.0 




0.1 
0.5 
1.0 

Field 


0.37l°;» 
0.881°;^? 

0.961°:?^ 


-8.691°:°? 

-8.46l°;°| 
-8.681°;°^ 
-8.47l°:°l 


, ,,7+0.045 
l-^l'-0.045 

1 lQO+0.011 

l-l-'°-0.()ll 
1 209+0-021 

1 107+0.001 


3 
11 
15 
105 



6. PROPERTIES OF GALAXIES AT FIXED TYPE, FIXED 
ENVIRONMENT, AND FIXED STELLAR MASS 

Given the covariance between galaxy stellar mass and 
environment, it is clear that in order to ascertain which 
has the largest causal effect on galaxy evolution, the 
properties of galaxies need to be measured func- 
tion of environment at fixed stellar mass, and vice versa. 

In addition to this, in § 5.2 we showed that the isF was 
also a function of stellar mass and environment. This 
means that in addition to any correlations induced from 
stcllar-mass-cnvironmcnt covariance, it is possible that 
part (or all) of the trends of the mean properties of galax- 
ies with environmental density (see Figures 6 and 7) may 
be caused only by the changing fsF, and not by a change 
in the properties of the star-forming galaxies and quies- 
cent galaxies themselves. This is an important distinc- 
tion, because if stellar mass or environment are causally 
linked to the quenching of star formation, we would ex- 
pect to measure a change in the SSFRs of star-forming 
galaxies with those parameters, beyond that seen simply 
from the changing isF with environment. 

Therefore, in order to untangle the interdependency 
of stellar mass, environment and isF, in this section we 
first examine how the isF depends on environment at 
fixed stellar mass, and vice versa. Thereafter we sep- 
arate star-forming and quiescent galaxies and measure 
how their SSFRs and D„(4000)s vary with environment 
at fixed stellar mass, and vice versa, in order to see which 
of these leaves the largest imprint on their stellar popu- 
lations. 

6.1. The Effect of Environment and Stellar Mass on the 
Star-Forming Fraction 

In the left panel of Figure 9 we plot the fsF as a func- 
tion of stellar mass in four environments: The field, R 
> 700 kpc, 200 < R < 700 kpc, R < 200 kpc. In the 



TABLE 5 
Galaxy Parameters for Quiescent 
Galaxies as a Function of Mass and 
Environment 



Radius D„(4000) Nstack 

Mpc 

(1) (2) (3) 



Log(M,/M0) > 10.7 



0.1 
0.5 
1.0 
Field 


1 794+0. 003 
1- '^^-0.003 
1 77Q+0.005 
1- ' "'-0.005 
1 70C+O.OO6 

1 734+0012 

^■'•5''_0.012 


57 

50 
27 
25 


10.0 < Log(M,/M0) < 10.7 


0.1 
0.5 
1.0 
Field 


1 r;7c: + 0.009 
1 p;Qa+0.003 
1 fiOQ+OOlO 

1 C7Q+O.O2I 
l.a(»_g 


25 
74 
15 
18 


9.3 < Log(M«/M0) < 10.0 


0.1 
0.5 

1.0 
Field 


1 4sq+0-037 
1 44S+0-011 

1.44CS_Q gjj 

1 QOC+0.004 

1.0Z0_Q 

1 9,.+0.064 
l-^'-'^-0.064 


3 
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3 
4 



right panel we plot the fsF as a function of clustercentric 
radius in three bins of stellar mass: LogM^/M© > 10.7, 
10.0 < LogM*/M0 < 10.7, and 9.3 < LogM^/Mo < 
10.0. Those arc slightly larger mass and environmental 
bins than those used in § 5, but given the additional 
separation of star- forming and quiescent galaxies, these 
adjustments are necessary to ensure adequate numbers 
of both in all stellar mass and environmental bins. 

Figure 9 shows that even at fixed environment, the fsF 
is still correlated with stellar mass, with more massive 
galaxies being less likely to be star forming. Interest- 
ingly, the relative decline in the fsF as a function of 
increasing stellar mass is consistent with being the same 
in all environments, roughly a factor of ~ 2-3 in all four 
environments. 

Although the relative decline in isF with increasing 
stellar mass occurs independently of environment, the 
absolute number of star-forming galaxies does depend on 
environment and this is demonstrated in the right panel 
of Figure 9 which shows the fsF as a function of environ- 
ment in the three stellar mass bins. The fsF is correlated 
with environment for galaxies over a wide range in stel- 
lar mass. The relative decline in fsF from the field to 
the cluster core is a factor of ~ 2-4, and within the un- 
certainties is consistent with being independent of stellar 
mass. The trend in decreasing tsF with increasing en- 
vironmental density appears to be slightly stronger than 
the trend with increasing stellar mass, suggesting envi- 
ronment may control this ratio slightly more than stellar 
mass; however, the uncertainties are still too large to de- 
termine if one is more dominant than the other. 

It is then quite clear from this analysis that both envi- 
ronment and stellar mass independently play a causal role 
in determining the fraction of star-forming galaxies at z 
~ i. It also appears that the quenching effects of stellar 
mass (self-quenching) and environment (environmental- 



17 



c 
.9 

'■4—' 

O 
CO 

LL 
LL 
C/) 




10^° 10^^ 0.0 0.5 1.0 1.5 

Stellar Mass (MJMJ Radius (Mpc) 

Fig. 9. — Left Panel: The fraction of star-forming galaxies (fsp) a function of galaxy stellar mass for galaxies in different environments. 
Right Panel: The igp as a function of environment for galaxies with different stellar masses. The fractional decline in Igp as a function of 
environment is a factor of ^ 3 at all stellar masses. Conversely, the fractional decline in isF as a function of stellar mass is a factor of ~ 2 
in all environments. This shows that the effects of stellar mass and environment on the fgp are independent and separable at z ~ 1. 
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Fig. 10. — Left Panel: Logarithm of the SSFR of star-forming galaxies as a function of stellar mass for galaxies in different environments. 
Right Panel: Logarithm of the SSFR of star-forming galaxies as a function of environment for galaxies with different stellar masses. The 
SSFR is correlated with stellar mass in all environments; however, but it is independent of environment at all stellar masses. This suggests 
that the primary factor in determining the SFRs of star-forming galaxies is their stellar mass, not their environment. 

lar mass, and vice versa. Our data confirm that this 
separabihty of environment and stellar mass on the 
continues to hold quite clearly up to z ^ 1. 



quenching) on the fsF are separable. Whatever process 
by which galaxies self-quench at a given stellar mass ap- 
pears to be equally efficient regardless of which environ- 
ment those galaxies are found in. Likewise, whatever 
process by which environment quenches star formation in 
galaxies appears equally efficient regardless of the stellar 
mass of the system. 

The id entical con c lusion was reached in the recent 
study bv lPeng et al.l (|2010( ) who showed using both the 
SDSS at z ~ 0, and zCOSMOS at z ~ 0.6 that the rel- 
ative increase in the quenched fraction of galaxies with 
increasing environmental density is independent of stel- 



6.2. The Specific Star Formation Rates of Star-forming 
Galaxies at Fixed Environment and Fixed Stellar 
Mass 

Now that we have established that both environment 
and stellar mass play an important and causal role in the 
quenching of galaxies at z ~ 1, we would like to under- 
stand the details of how those processes might work. In 
this subsection we consider the properties of galaxies of 
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Fig. 11. — Left Panels: D„(4000) of quiescent and star-forming galaxies as a function of stellar mass for galaxies in different environments. 
Right Panels: D„(4000) of quiescent and star-forming galaxies as a function of environment for galaxies with different stellar masses. The 
Dn(4000) of both quiescent and star-forming galaxies is a function of stellar mass in all environments; however, it is independent of 
environment at all stellar masses. Given the correlation between D„(4000) and the age of the stellar population, this suggests that stellar 
mass is the primary factor in determining the age and/or metallicity of a galaxy, not its environment. 

each type individually - either star-forming or quiescent 
- at fixed environment, and at fixed stellar mass. As in 
§ 5, this analysis is based on measurements made from 
stacking galaxies within different stellar mass and envi- 
ronmental bins. 

In the left panel of Figure 10 we plot the SSFRs of the 
star-forming galaxies in different environments as a func- 
tion of stellar mass and in the right panel we plot their 
SSFRs as a function of clustercentric radius in different 
stellar mass bins. Figure 10 shows that the SSFR of star- 
forming galaxies is correlated with their stellar mass, but 
this correlation appears to be basically independent of 
the environment they live in. Given the clear dependence 
of the isF on both stellar mass and environment, this is 
a surprising result. If both stellar mass and environment 
are responsible for quenching star formation it might be 
expected that the SSFRs of star-forming galaxies should 
decrease both as a function of environment and stellar 
mass, not just stellar mass alone. 

A lack of dependence of the SSFR of star-forming 
galaxies on environment is also observed in the local uni- 



verse. iBalogh et ahl ()2004aD found that in the 2dFGRS 
and SDSS the galaxy population is bi-modal in terms of 
the EW(Ha). They found that galaxies with EW(Ha) < 
4A dominate the high density regions, and galaxies with 
EW(IIa) > 4A dominate the low density regions. Sim- 
ilar to our results at z ^ 1, the fraction varies strongly 
with environment, but the mean EW(Ha), which is di- 
rectly related to the SSFR, is n early independent of en - 
vironment. Independently, both lKauffmann et al.l ()2004l ) 
and Peng et al. (201Q|) have found the same effect in the 
SDSS. They use different definitions of star-forming and 
quiescent galaxies, as well as different measures of the 
SSFR, which suggests the result does not strongly de- 
pend on those choices. 

The most likely, and perhaps only reasonable expla- 
nation for a lack of a dependence of the SSFR of star- 
forming galaxies on environment is that the timescale 
over which environment quenches star formation in 
galaxies is very rapid. If so, environmental quenching 
will move galaxies out of the star-forming classification 
and into the quiescent classification before a drop in their 
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SSFRs is measured. Such a process may have a measur- 
able effect on the D„(4000) of quiescent galaxies as a 
function of environment as well as the number of post- 
starburst galaxies as function of environment, issues we 
investigate further in § 6.3 and § 7. 

6.3. The Dn(4000) of Star- Forming and Quiescent 
Galaxies at Fixed Environment and Fixed Stellar 
Mass 

In Figure 11 we plot the D„(4000) measures from the 
stacked spectra as a function of stellar mass at fixed en- 
vironment (left panels) and as a function of environment 
at fixed stellar mass (right panels) . We can now consider 
the properties of both the star-forming and quiescent 
galaxies, which are shown in the bottom and top pan- 
els of Figure 11, respectively. We have also listed these 
parameters in Tables 4 and 5. The trends in D„(4000) 
with stellar mass and environment are qualitatively sim- 
ilar to the trends in the SSFR with these parameters. 
For both the star-forming and quiescent galaxies the 
D„(4000) measures are strongly correlated with stellar 
mass, regardless of environment, and the D„(4000) ap- 
pears to be independent of environment, regardless of 
the stellar mass. The exception is the lowest-mass qui- 
escent galaxies, which suggest a trend in D„(4000) with 
environment. This trend is quite interesting if real; how- 
ever, low- mass quiescent galaxies are rare at z ~ 1, and 
the statistics in these stacks are still somewhat limited. 
Even though the bootstrap error bars are small, those 
bins contain between 3-15 galaxies, and more data are 
probably required to confirm this interesting potential 
trend. 

Together, Figures 10 and 11 comprise one of the main 
conclusions of this paper, namely that it appears that 
at z ~ 1 stellar mass is the parameter that is primar- 
ily responsible for determining the stellar populations of 
both star-forming and quiescent galaxies, not their envi- 
ronment. Instead, the main role of environment is that 
it determines the fsp, and also at some level, the mean 
stellar mass of galaxies. Once the covariance of isp with 
environment is controlled for, as well as the covariance of 
stellar mass with environment, it appears that environ- 
ment does not affect the mean properties of both star- 
forming and quiescent galaxies, only the relative fraction 
of each. 

In § 6.2 we suggested that the lack of a dependence 
of the SSFR of star-forming galaxies on environment 
could be explained by a rapid environmental quenching 
timescale that quickly transforms star-forming galaxies 
into quiescent galaxies. If so, we would also expect that 
the D„(4000) of the star- forming galaxies would be in- 
dependent of environment, as is observed. However, we 
might also expect that the D„(4000) of the quiescent 
galaxies would correlate with environment, because a 
higher fraction of galaxies in high-density regions have 
been quenched by their environment, hence they may 
have different D„(4000) at fixed stellar mass compared 
to lower-density environments. Except for the lowest- 
mass bin, where the data are fairly sparse, we do not 
see evidence for such a dependence. As we will show 
with some simple quenching models in § 8, the lack of 
a correlation between the D„(4000) of quiescent galax- 
ies and their environment is still possible under certain 
circumstances, even if there is a substantial amount of 



environmental quenching. 

7. THE POSTSTARBURST POPULATION 

Our current hypothesis is that the timescale for envi- 
ronmental quenching is rapid enough to transform star- 
forming galaxies into quiescent galaxies before a decline 
in their SSFRs can be observed. One way to test this 
is to try to identify the population of recently-quenched 
galaxies and see how it correlates with both environment 
and stellar mass. Defining a complete sample of recently- 
quenched, or "poststarburst" galaxies is difficult because 
the quantitative description of this class of galaxies has 
varied within the literature. Typically, when rest-frame 
optical spectroscopy is available, many authors have 
considered galaxies classifie d as "K-|- A galaxies" (e.g^ 
Balogh et al. 1999; Dressier et al.|[T99 9: Pqggianti ct aT] 
1999t' ILe Borgne et all 120061: iPoggianti'et'al.i i200l 
Yan et al.ll2009() as poststarburst galaxies or regular star- 
forming galaxies whose star formation was abruptly ter- 
minated before the epoch of observation. The quanti- 
tative definition of K-l-A in itself has taken a range of 
values in different studies; however, in general it is the 
subset of galaxies with weak [Oil] emission and strong 
US absorption-line strength. The lack of [Oil] emission 
is an indicator of a lack of ongoing star formation, and 
the strong US absorption is indicative of a population 
of A-stars which have lifetimes of ~ 1 Gyr. Therefore, 
K-t-A galaxies can be considered to have no ongoing star 
formation but a young stellar population less than 1 Gyr 
old. 

The GCLASS spectra are of sufficient quality that the 
EW([OII]) can be measured in individual spectra as weak 
as l-SA (see § 5); however, KS is a weak absorption fea- 
ture that at z 1 is in a wavelength range contaminated 
with numerous sky lines, making it difficult to use the 
common K-l-A definition with these data. Fortunately, 
at z 1 we can take advantage of the relatively young 
absolute age of most galaxies and the rapid non-linear 
evolution of D„(4000) as a function of galaxy age. As 
Figure 4 shows, the D„(4000) of a single-burst popula- 
tion evolves from a value of D„ (4000) — 1.0 at zero-age to 
D„(4000) = 1.5 at one Gyr after the initial burst. There- 
after the strength of the break begins to evolve much 
more slowly with time. In principle, a weak D„(4000) 
can be used as a proxy for a young stellar population, 
much in the same way that the K-|~A definition uses US 
a s a proxy for a young stellar population. 

iBalogh et aO Il999h showed using synthetic spectra 
tracks that galaxies that have had their star formation 
abruptly truncated show a correlation between EW(H(5) 
and D„(4000) (their Figure 11). They showed that for 
galaxies with EW(II5) > 5A (the typical definition for a 
K-t-A galaxy) the typical range of D„(4000) values is 1.0 
< D„(4000) < 1.45. Therefore, it might be expected that 
selecting galaxies with, 1) no detectable [Oil] emission, 
and 2) 1.0 < D„(4000) < 1.45, these may be analogous 
to the standard K-l-A selection with the caveat that this 
selection will miss older galaxies that experienced a re- 
juvenation period that was subsequently truncated. 

We tested what type of stellar population this selec- 
tion criteria selects by stacking the spectra of all galaxies 
that met the criteria. This stacked spectrum is plotted in 
Figure 12. The spectrum is dominated by Balmer lines 
and is clearly similar to the spectra oiK+A galaxies. The 
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Fig. 12. — Mean stacked spectrum of candidate poststarburst galaxies with 1.0 < D„(4000) < 1.45 and no detectable [Oil] emission. The 
spectrum has been smoothed with a SA boxcar. Prominent absorption (emission) features are marked in red(green). The stacked spectrum 
clearly shows the full Balmer series of absorption lines. It has an EW(H(5) = 5.2 ± 0.6, and EW([OII]) = 0.9 ± 0.7 and is classified as 
a K+A galaxy. The stacked spectrum suggests that on average the [Oil]- and D„ (4000)-based poststarburst selection technique selects 
similar galaxies as the more typical K+A selection. 
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Fig. 13. — Left Panel: Fraction of poststarburst galaxies as a function of stellar mass in different environments. Right Panel: Fraction of 
poststarburst galaxies as a function of environment in different stellar mass bins. The number of poststarbursts increases with decreasing 
clustercentric radius, except for the most massive galaxies, of which few are poststarbursts in any environment. Likewise, pos tstarbursts 
appea r to be rare (< 2% of the population) at all stellar masses in the field comparison sample. The DEEP2 K+A selection from lYan et al.l 
II2009I ) is shown for comparison and agrees well with our result. 

now plot the fraction of poststarburst galaxies as a func- 
tion of stellar mass in different environments in the left 
panel of Figure 13. Likewise, in the right panel we plot 
the fraction of poststarburst galaxies as a function of en- 
vironment in several bins of stellar mass. As a reference 
we denote the DEEP2 measurement of the K+A fraction 
at z ~ 0.8 from Yan et al., (.2009,) as the red point. This 
value is consistent with our field measurement, and fur- 
ther support the claim that our poststarburst selection 
criteria selects similar galaxies as the K+A criteria. 

Figure 13 shows that the fraction of poststarburst 
galaxies is simultaneously a function of stellar mass, and 
environment. Less massive galaxies are more frequently 
poststarbursts compared to massive galaxies in all envi- 
ronments, and galaxies in the highest-density environ- 
ments such as the cluster cores, are more frequently 
poststarbursts. The latter statement is not true for the 
LogAf0 > 10.7 galaxies, which have almost a zero post- 
starburst fraction in the clusters. 

It is interesting that poststarbursts are more common 
in the high-density cluster environment at z ^ 1 com- 



stack has an EW(H5) = 5 2A ± 0.6A, and an EW([OII]) 
= 0.9A ± 0.7A. Using the lBalogh eTaLl (|1999f ) definition 
of K+A galaxies: EW(H(5) > 5A, and EW([OII]) < SA, 
the mean galaxy selected with the D„(4000)-based and 
EW([OII])-based criteria is in fact a K+A galaxy, sug- 
gesting this selection is also a good criteria for a post- 
starburst galaxy. We note that this criteria almost cer- 
tainly does not a select the complet e sample of post - 
starbursts. As has b een discuss ed bvlYan et al.l ()2009[ ). 
iLemaux et all (j2010l ) and Koc evski et al.l (|2010f ). even 
the more rigorous K+A definition misses poststarburst 
galaxies with [Oil] emission from an AGN, and can also 
misclassify star-forming galaxies as poststarbursts when 
the [Oil] line is highly extinguished by dust. In our anal- 
ysis we are primarily interested in the environmental and 
stellar-mass dependence of the poststarburst population. 
Even if the definition selects only a subset of all poststar- 
bursts, provided it does so consistently the comparison 
of the poststarburst fraction at different stellar masses 
and environments should still be meaningful. 
With a clear definition of a poststarburst galaxy, we 
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pared to the field. Taking an average over galaxies with 
9.3 < LogM,/Mo < 10.7 shows that they are 3.1 ± 1.1 
times more common in galaxy clusters compared to the 
field. Not only is this excess of poststarburst galaxies 
in the high-density cluster environment more direct evi- 
dence for the importance of environmental quenching of 
star formation at z ~ 1, it also shows that environmental 
quenching is sufficiently abrupt to create the poststarburst 
signature in galaxies. 

A similar overdensity of poststarburst galaxies in and 
aroun d clusters has also been seen by iPoggianti et al.l 
(|2009[) in the EDisCS cluster sample at 0.4 < z < 0.8. 
They find that the K-|-A fraction is a factor of ~ 2 higher 
in cluster s compared to a con trol sample of field galaxies. 
Likewise. Balo gh et al.l ()2011l) have also seen evidence for 
a large population of "green valley" galaxies in galaxy 
groups at z ~ 1. Taken together, these results all consis- 
tently point to high-density regions playing an important 
role in quenching galaxies up to 2 1. 

Indeed, given that self quenching should occur in all 
environments, if we take Figure 13 at face value it seems 
to imply that the primary cause of the poststarburst 
phase is environmental quenching. Quantitatively the 
two match as well, as the both the fraction of quiescent 
galaxies and the poststarburst fraction increase by a fac- 
tor of ~ 3 as we move from the field to the cluster environ- 
ment. We have speculated that a rapid environmental- 
quenching timescale can explain the lack of a correlation 
between the SSFRs of star-forming galaxies and their en- 
vironment, and the increasing fraction of poststarbursts 
with increasing environmental density lends further sup- 
port to this hypothesis. Our results also suggest that the 
environmental-quenching timescale may be more rapid 
than the stellar-mass quenching timescale. The correla- 
tion between the SSFR of star-forming galaxies and their 
stellar mass (§ 6.2) suggests that at some level, galaxies 
self-regulate their growth by winding down their SSFRs 
as they grow in stellar mass. Environment appears to be 
a stochastic process that can quickly truncate the star 
formation in a galaxy once it has fallen into a sufficiently 
large overdensity. 

Another interesting point from Figure 13 is that al- 
though poststarbursts are substantially more common in 
the cluster environment, they still make up a relatively 
small fraction of the overall population, even for lower- 
mass galaxies. Galaxies that fit our poststarburst defi- 
nition are 10% ± 2% of all cluster galaxies with 10.7 < 
LogMH./MQ < 9.3, and 0.5% ± 0.5% of cluster galaxies 
with LogM> 11.0. If we compare the number of post- 
starbursts to the number of quiescent galaxies in the clus- 
ter, and assume that environment is the primary source 
of the poststarburst population, it suggests that 16% ± 
3% of quiescent galaxies with 9.3 < LogM^/Mo < 10.7 
have had their star formation truncated by their envi- 
ronment within the last ^ 1 Gyr. It also suggests that 
very few of the quiescent galaxies with LogM*/M0 > 
10.7 have been truncated by their environment within 
the last Gyr. This recently-quenched fraction is mod- 
est and suggests that even if environmental quenching 
dominates the growth of the low-mass end of the quies- 
cent population, the quenching rate for such galaxies is 
still fairly low at z 1, and that the vast majority of 
quenched galaxies at that epoch were quenched at least 



> 1 Gyr prior (i.e., z > 1.3). 

8. A SIMPLE QUENCHING MODEL TO EXPLAIN THE 
LACK OF A D„{4000)-ENVIRONMENT CORRELATION 
AT FIXED STELLAR MASS 

The existence of a poststarburst population that is 
associated primarily with high-density regions appears 
to support a rapid environmental-quenching timescale, 
which in turn provides a good explanation for the lack 
of a correlation between the SSFRs and D„(4000)s of 
star-forming galaxies and their environment. However, 
in § 6.3 we also showed that the D„(4000) of the quies- 
cent galaxies is independent of their environment. We 
know from both the strong decline in the isp with in- 
creasing environmental density, as well as the increase in 
the poststarburst fraction with increasing density that 
environmental processes must be responsible for quench- 
ing a significant fraction of the quiescent population. It 
seems then difficult to understand why the D„(4000) of 
quiescent galaxies does not correlate with environment 
at fixed mass, but instead correlates with stellar mass 
at fixed environment. In this section we construct some 
simple quenching models to investigate how much the 
D„(4000) of quiescent galaxies can differ between envi- 
ronments in different quenching scenarios. 

8.1. A Basic 0^(4000) Evolution Model 

On purpose, our model is simple and designed to pro- 
vide a basic illustration of the problem, not a detailed 
quantitative prediction for the evolution of D„(4000) in 
different environments. Our approach is to simulate a 
population of star-forming galaxies, assuming they are all 
the same stellar mass, and then quench them at different 
rates to follow the evolution of the average Dn(4000) of 
only the quenched galaxies at fixed stellar mass. To make 
this model requires that we have a prediction of how the 
D„(4000) evolves with time for both star- forming and 
quiescent galaxies. 

We assume that star-forming galaxies form stars at a 
constant rate until they are quenched, either by their 
environment, or by some self-regulating process related 
to their stellar mass. As the star-formation histories 
in Figure 3 show, galaxies with a constant SFR have 
a slowly evolving D„(4000) that remains near a value 
of ~ 1.2. Therefore, for simplicity, we assume that the 
D„(4000) of star-forming galaxies is exactly equal to 1.2 
at all times while they are star forming. In terms of 
the quiescent galaxies, the red track in Figure 3 shows 
that the D„(4000) of a single-burst population evolves 
rapidly after the initial burst, and thereafter the rate of 
evolution begins to slow. We will assume that once star- 
forming galaxies have been quenched, they follow along 
the D„(4000) track for a passively-evolving single-burst 
galaxy. We note this is not a completely fair represen- 
tation, as the D„(4000)'s of galaxies that have been star 
forming for longer periods will evolve slightly slower than 
single-bursts; however, given that we are more interested 
in illustrating the difference in D„(4000) between differ- 
ent environments and different quenching scenarios, not 
predicting its precise value, a more complex treatment of 
the star formation histories is not warranted. 

As a first step in the model, we define three scenarios 
by which the quenching rate of galaxies might evolve, 
regardless of whether environmental quenching or self 
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Time (Gyr) Time (Gyr) 

Fig. 14. — Left Panel: Model tracks of the fraction of quenched galaxies as a function of time for simple analytic functions for the quenching 
rate; linear (green), exponentially declining (red), exponentially increasing (blue). Right Panel; Evolution of the average D„(4000) for the 
quiescent galaxy population using the three quenching rate functions. A galaxy that is quenched at t = and hence is maximally old is 
shown for reference. The exponentially declining quenching rate has the strongest average Dn(4000) because more galaxies quench early 
and hence have more time to evolve passively. An exponentially increasing quenching rate has the weakest average Dn(4000) because the 
population of quiescent galaxies is dominated by galaxies that were recently quenched. 




Fig. 15. — Left Panel; Model tracks for the fraction of quenched galaxies as a function of time using the analytic quenching functions in 
Figure 14. These models include separate mass-quenching components (dotted lines), and a environmental-quenching components (solid 
lines). The red and blue models have the same amplitude for mass and environmental quenching but the quenching rates evolve differently 
with time. The orange and green models have different amplitudes for the mass and environmental quenching components, but the rates 
evolve the identically with time. Right Panel; Relative difference in the inferred D„(4000) for the model that has both environmental and 
mass quenching as in the left panel (simulating a cluster) compared to a model that has just mass-quenching (simulating the field). If the 
rate of mass and environmental quenching evolves identically, there is no difference in the D„(4000) between the cluster and field, even 
if the amplitudes of each are different. The model where environmental quenching dominates at early times (red) has a larger Dn(4000) 
in the cluster compared to the field. When environmental quenching dominates at late times (blue), the quiescent galaxies in the cluster 
have weaker Dn(4000) than the field, because they are dominated by recently-quenched galaxies. The arrows show the effect of increasing 
(purple) or decreasing (black) the importance of environmental quenching relative to mass quenching. 
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quenching is responsible. We consider the evolution of 
the quenching rate over a period of 6 Gyr, which is the 
age of the universe at z ~ 1. This is the maximum 
amount of time over which galaxies have had to quench, 
and will illustrate the maximum difference between the 
properties of galaxies in different environments that can 
be achieved. 

The first quenching rate we consider is where the 
quenching rate of galaxies proceeds linearly, and the rate 
does not evolve with time. In the left panel of Figure 14 
we plot the fraction of quenched galaxies as a function of 
time for this model in green. We also consider a model 
where most of the quenching occurs at early times, where 
we parameterize the quenching rate as a decaying expo- 
nential (red). Lastly, we consider a model where most 
of the quenching occurs at late times, where the quench- 
ing rate is parameterized as an increasing exponential 
(blue). The models are normalized so that 100% of the 
quiescent galaxies are quenched by t = 6 Gyr. We have 
chosen an exponential to parameterize the "early" and 
"late" quenching scenarios; however, we note that any 
sharply increasing or decreasing function would be suffi- 
cient. 

In the right panel of Figure 14 we plot the predicted 
evolution of the mean D„(4000) of quiescent galaxies as 
a function of time for the three quenching functions. For 
reference we also plot the mean D„(4000) of a purely 
passively evolving galaxy to show the upper limit of 
D„(4000) as a function of time. As expected, the mean 
D„(4000) of quiescent galaxies depends on the epoch 
when they are quenched. The exponentially declining 
quenching model has the largest mean D„(4000) be- 
cause most of the galaxies have quenched at early times 
and hence have had more time to age. The exponen- 
tially increasing quenching model has the smallest mean 
D„(4000), because at all times the quiescent population 
is dominated by galaxies that have only recently been 
quenched. 

8.2. A Model for the Competing Effects of 
Environmental and Self Quenching on Dn(4000) 

We can now add complexity to the model by assuming 
two independent quenching mechanisms: one related to 
stellar mass, which we will call self quenching, the other 
related to environment, which we refer to as environ- 
mental quenching. In particular, we would like to under- 
stand the difference in the D„(4000) of quiescent galaxies 
that results from a model that has both self quenching 
and environmental quenching - which should represent 
a high-density region such as a cluster - as compared to 
the predicted D„ (4000) of quiescent galaxies from models 
that have only self quenching - which should represent 
the low-density regions found in the field population. 

There are a large range of quenching models that can 
be built to make this comparison; however, the most 
relevant things to test in terms of a difference in the 
D„(4000) of quiescent galaxies with environment is 1) 
the amplitude of self quenching as compared to environ- 
mental quenching, i.e., which process is most important 
for creating the quiescent population, and 2) how the 
time evolution of the self-quenching and environmental- 
quenching rates affects the D„(4000)s. 

In the left panel of Figure 15 we plot the fraction of 
quenched galaxies as a function of time for models that 



include both self quenching and environmental quench- 
ing. Each model for the quenching rate is plotted in a 
different color, and within those colors the fraction of 
galaxies quenched by environment is shown as the solid 
curve and the fraction of galaxies that are self quenched 
is shown as the dotted curve. 

We test four models to explore the importance of items 
1) and 2) above. Firstly, we examine the case where 
environmental quenching and self quenching have equal 
importance, each quenching 50% of the population, but 
the rate at which they quench galaxies evolves differ- 
ently with time. The red model shows the case where 
environmental quenching dominates at early times, and 
self quenching dominates at late times. The blue model 
shows the converse, with self quenching dominating at 
early times and environmental quenching dominating at 
late times. 

In the right panel of Figure 15 we plot the difference 
in the evolution of the mean D„(4000) of the quiescent 
galaxies for these models that include both self and en- 
vironmental quenching compared to the mean D„(4000) 
for a model that just has self quenching (with the same 
self-quenching rate evolution as in the combined model) . 
The red line shows that when environmental quenching 
dominates early, quiescent galaxies in high-density envi- 
ronments will have stronger D„(4000)s than those in low- 
density environments. This occurs because the galaxies 
in high-density environments quench earlier on average 
and hence have more time to evolve. The blue curve 
shows that when environmental quenching dominates at 
late times, quiescent galaxies in high-density regions ac- 
tually have weaker D„(4000) than those in low-density 
regions. This occurs because high-density regions have a 
larger population of recently-quenched galaxies that are 
young and move the mean of D„(4000)s to a lower value. 

The blue curve is interestin g and is worth co mment- 
ing on. A prediction from the lPeng et al.l (j2010( ) empir- 
ical quenching model is that self quenching should dom- 
inate the early evolution of galaxies and environmental 
quenching should dominate their later evolution as the 
growth of structure in the universe proceeds. Our sim- 
ple model for the D„ (4000) evolution of quiescent galax- 
ies shows that if this is the case, we might expected to 
see something quite counter-intuitive, that at fixed ste- 
lar mass, quiescent galaxies in clusters should actually 
be younger than t hose in the fie ld. In general, th e con- 
verse is seen jejr., Tho mas et allf 2005: Nclan ct al.' 2001 
Gallazzi et al.ll2005: van Dokkum & : van jlexMarelJMS 
Thomas et al.ll2010[ ). which would seem to argue against 
environmental quenching dominating at late times. In- 
deed, if anything, it would seem to suggest that envi- 
ronmental quenching should dominate at earlier epochs. 
On the other hand, we should also consider that most of 
those studies have focused on primarily elliptical galax- 
ies, not all quenched galaxies. In general clusters contain 
a much higher fraction of SO galaxies than the field (e.g., 
iDressleiillOS O") and it may be that when SO's are included 
with the elliptical galaxies, the population of quenched 
galaxies in clusters (E -I- SO) is in fact younger than the 
population of quenched galaxies in the field (dominated 
by E). 

Although it is unclear how to reconcile this issue, it 
is important to point out that the absolute difference in 
the D„(4000)s between our simulated cluster and field is 
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quite small, with AD„(4000) = 0.1 to 0.2. Our model 
assumes an equal contribution from self quenching and 
environmental quenching. In principle, the AD„(4000) 
between cluster and field can be increased or decreased 
by adjusting the amplitude of self quenching compared 
to environmental quenching. Increasing the importance 
of self quenching over environmental quenching decreases 
the difference between cluster and field, and likewise, in- 
creasing the importance of environmental quenching over 
self quenching increases the difference between cluster 
and field. This shows that a lack of a difference in the 
D„(4000) of quiescent galaxies as a function of environ- 
ment would be expected if self quenching dominates over 
environmental quenching at all epochs. Indeed, if greater 
than two-thirds of the cluster galaxies are self quenched, 
this would produce a AD„(4000) that is < 0.1, and hence 
is unlikely to be detected with our data. 

While it is the simplest answer, a dominance of self 
quenching over environmental quenching is not a fully 
satisfactory explanation for the lack of a D„(4000)- 
environment correlation. We know that the change in 
the isF with environment is similar to the change in the 
isF with stellar mass, and that there are a substantial 
number of poststarbursts associated with clusters. Both 
of these suggest that environmental quenching is at least 
as important in the transformation of galaxies as is self 
quenching. In order to explore this further, we now con- 
sider a second set of models. This time where the am- 
plitude of self qucnciiing and environmental quenching 
are quite different, but where the quenching rates evolve 
in a similar way. For this comparison we use the linear 
quenching rate, where in fact the rate does not evolve 
with time, but we note that the exponential models that 
do evolve with time provide a similar comparison. 

The green curves in Figure 15 show the model where 
environmental quenching dominates, quenching 80% of 
the population by t = 6 Gyr, and the orange curve 
shows the converse, where self quenching dominates and 
quenches 80% of the population by t = 6 Gyr. Interest- 
ingly, the AD„ (4000) between cluster and field for these 
models is precisely zero. This demonstrates that it is 
not critical that self quenching dominates over environ- 
mental quenching in order to produce similar D„(4000) 
in high-density and low-density regions. All that is re- 
quired is that the self quenching rate evolves in the same 
way as the environmental quenching rate. If they do, 
then no matter the amplitude of each, the AD„(4000) 
between cluster and field is precisely zero. Put another 
way, the D„(4000) depends on when galaxies quench 
in cosmic time, not what is the source of the quench- 
ing. Provided that the environmental-quenching and 
self-quenching rates evolve so as to quench galaxies simi- 
larly at all epochs, the fraction of quenched galaxies will 
be higher in clusters; however, then the mean D„(4000) 
of the quiescent galaxies will be the same in all environ- 
ments. 

Putting both sets of models together provides two 
plausible explanations for why the D„(4000) of quies- 
cent galaxies does not depend on environment. First, 
if self quenching dominates over environmental quench- 
ing, then the D„(4000) will be the same in all envi- 
ronments, regardless of how the two quenching rates 
evolve. Alternatively, provided that the self-quenching 
and environmental-quenching rates evolve similarly with 



time, then regardless of which process dominates no dif- 
ference in the D„(4000) of quiescent galaxies with en- 
vironment will be seen. The latter would be a better 
explanation for our data, which suggests that environ- 
mental quenching is roughly equally important as mass 
quenching. 

9. SUMMARY 

In this paper we have studied the correlations between 
the properties of galaxies and their environment and stel- 
lar mass. We first considered the average properties of 
a mass-limited sample of galaxies with LogM*/M© > 
9.3 (§ 5), and then subdivided the galaxies into differ- 
ent types (star-forming or quiescent) and studied the 
properties as a function of environment at fixed stellar 
mass, and vice versa. We also identified the population of 
recently-quenched galaxies and studied their abundance 
as a function of environment and stellar mass. Lastly, we 
constructed several simple quenching models to explain 
some of the (lack of) correlations in the data. The main 
results of the paper are, 

1. For galaxies with LogM^/M© > 9.3, the well-known 
correlations between fs^, SSFR, SFR, D„(4000), and 
color with environment seen in the local universe are al- 
ready clearly in place by 2; ~ 1. 

2. For cluster galaxies, the well-known correlations be- 
tween properties such asisF, SSFR, SFR, D„(4000), and 
color with stellar mass seen in the local universe are also 
in place at 2; ~ 1. 

3. Similar to the local universe, there is a covariance 

between the mean stellar mass of galaxies and their en- 
vironment, with more massive galaxies living in higher- 
density environments. 

4. The relative decline in the isp with environment at 
fixed stellar mass is a factor of ~ 2-4 and the relative 

decline in the tsF with stellar mass at fixed environment 
is a factor of ~ 2-3. Both stellar mass and environment 
affect the fsF, and their relative effects are separable up 
to 2; ~ 1. 

5. Controlling for the covariance in stellar mass and fsp 
with environment we find that the SSFRs of star-forming 
galaxies are correlated with stellar mass at fixed environ- 
ment, but is nearly independent of environment at fixed 
stellar mass. 

6. The D„(4000)s of both star-forming and quiescent 

galaxies are correlated with their stellar mass at fixed 
environment, but is nearly independent of environment 
at fixed stellar mass. 

7. The data suggest that stellar mass is the primary pa- 
rameter that determines the stellar populations of galax- 
ies. Environment plays a different role, controlling the 
fsF, and to a lesser degree, the mean stellar mass of 
galaxies. 

8. The population of poststarburst galaxies is correlated 
with both environment and stellar mass, in that post- 
starbursts are more common amongst lower-mass galax- 
ies and in higher-density environments. Poststarburst 
galaxies with 9.3 < LogM^/M© < 10.7 are 3.1 ± 1.1 
times more common in clusters at z ~ 1 compared to 
the field. 

9. The lack of a correlation of the SSFRs and D„(4000) 
of star-forming galaxies and their environment is best- 
explained by a rapid environmental quenching timescale. 
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The rapid environmental quenching tiniescale is further 
supported by the clear association of poststarbursts with 
high-density environments. 

10. Despite an overabundance of poststarbursts in clus- 
ters compared to the field, poststarbursts still make a 
relatively small fraction of the total cluster population. 
Only 10% ± 0.2% of cluster galaxies at z - 1 with 9.3 < 
LogM*/M0 < 11.0 are poststarbursts. 

11. Simple quenching models show that the lack of a cor- 
relation between the D„(4000) of quiescent galaxies and 
their environment can be explained two ways: Either self 
quenching dominates over environmental quenching at 
z > 1, or else the evolution of the self quenching rate mir- 
rors the evolution of the environmental quenching rate. 

Overall, this study of z ~ 1 clusters shows that 
the trends in galaxy properties with environment 
and stella r mass that are seen in the local uni- 
verse (e.g., 'Kauff mann et all 12004 iBalogh et al.ll2004bl: 
Baldrv et a l. 2006: iWeinmann et al.l l2006t iPeng et al.l 
2010| ) are already well-established at z ^ 1. It appears 
clear from our study as well as both lower- and higher- 
redshift studies that while the fraction of star-forming 
galaxies is a strong f unction of their environment, (e.g., 
Balogh et al]|2004bt iKauffmann et ahl 12004 iPeng et alj 
20ld iPatel et al.l 120111: lOuadri et al.l 120111) . the mean 
properties of the star-forming and quiescent galaxies at 
fixed stellar mass are not. This means that stellar mass 
is the primary parameter linked to determining the stel- 
lar populations of galaxies. 

Our simple modeling shows that the lack of corre- 
lation of the properties of star-forming and quiescent 
galaxies with their environment can be understood if 
the environmental-quenching timescale is rapid, and that 



the evolution of the self-quenching and environmental- 
quenching rates mirrors each other - regardless of which 
processes dominates the overall quenching process. This 
hypothesis is further supported by the substantial pop- 
ulation of poststarburst galaxies that are only found in 
the high-density cluster environment. 

The next obvious step to better understanding of the 
process of environmental quenching will be to connect the 
evolution of the cluster galaxy population over a larger 
redshift baseline. Doing so will allow us to differenti- 
ate between different quenching models and get a better 
handle on the physical processes behind environmental 
quenching. 
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